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EDITORIAL 


Com a publicação deste número duplo pretendeu-se e conseguiu-se reduzir o intervalo entre as 
sucessivas edições da Técnica. 


Contribuiu de modo decisivo para a recuperação agora iniciada, a imediata publicação dos 
originais cujas provas de revisão se deram por findas. Desta atitude, poderá resultar que a ordem de 
entrada dos artigos não seja exactamente a mesma que se observa na publicação, mas tal facto fica 
salvaguardado pelo anúncio da data de entrada do original. 


Se a afluência de artigos for numerosa e se distribuir harmonicamente pelas diversas secções e áreas 
cientificas, se a revisão das provas for rápida, o inconveniente apontado tenderá a desvanecer-se. 


Os aspectos focados pesam de modo importante nas limitações materiais que temos vindo a referir em 
números anteriores. Nesta matéria temos que ser intransigentes no que respeita à introdução de 
alterações no original inicial e já em trabalho de composição. 

Mais uma vez se solicita que se respeitem as NORMAS relativas à extensão dos artigos. 

No que se prende ao estritamente material, nomeadamente à forma e tempo de composição, 
continua-se a trabalhar. 


Relativamente à colaboração recebida e natureza de artigos considerados, a experiência colhida e uma 
prospecção no dominio de potenciais leitores aconselha-nos o seguinte: 


— Manter dois ARTIGOS CONVIDADOS por número. 
Nesta área temos recebido pronta e excelente colaboração. 


— Reservar na secção de ARTIGOS uma parte importante, pretenderiamos quase total, para os arti- 
gos de sintese que, em linguagem rigorosa mas acessivel, facultem aos estudantes as últimas 
novidades da Ciência e da Técnica, seu interesse actual e futuro. 


Cria-se assim um espaço novo que terá significativa projecção entre estudantes e até de outros 
que, embora especialistas, o não são nas diversas áreas contempladas. 


— Manter e incentivar as NOTAS CIENTÍFICAS. 


Trata-se de um espaço eminentemente cientifico vocacionado para apresentar resultados de 
investigação ou experiência colhida na realização de trabalhos técnico-cientiíficos. 


Aqui a colaboração tem rareado. 


RESUMINDO 


— Comos'“ ARTIGOS CONVIDADOS” a TÉCNICA reserva-se o direito de estimular o que pre- 
tende publicar; 


— Comos “ARTIGOS” espera uma larga comparticipação entre professores e alunos e de que por 
acções concertadas se construam originais de interesse; 


— Comas'“NOTAS CIENTÍFICAS” a TÉCNICA oferece um espaço de publicação à actividade 
de Investigação e Desenvolvimento do 1.5. T. e não só. Com esperança, continua aguardando. 


M. de Abreu Faro 
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SUMÁRIO 


O controlo da gestão de stocks de um armazém que abas- 
tece uma linha de montagem numa unidade industrial de 
produção em série é considerado através da interacção dis- 
tribuida entre os sistemas intervenientes. As actividades de- 
senvolvidas são identificadas fundamentalmente em dois 
grupos: comunicação, representando trocas de mensagens 
entre os sistemas componentes, e processamento, corres- 
pondendo a acções intensivas de tratamento da informa- 
ção. O artigo detalha o diagrama temporal, analisando as 
actividades locais e as mensagens de comunicação trocadas 
entre os diversos sistemas componentes. Processamentos 
típicos são referidos, apresentando-se algoritmos de selec- 
ção de alternativas nomeadamente um baseado na cha- 
mada regra de stock minimo. 

O artigo estuda a sistematização de actividades em am- 
biente que designaremos por não estruturado. A não estru- 
turação é por exemplo resultado da intervenção de opera- 
dores humanos. A definição do próprio sistema global, 
envolvendo a especificação de comportamentos complexos, 
é uma tarefa que subjaz o presente trabalho, mas que não é 
aqui considerada. Introduzindo progressivamente detalhes 
complementares, o artigo identificará no final aspectos es- 
senciais do problema. 

O trabalho exposto é parte integrante de um projecto 
operacional numa empresa industrial. A redução da inter- 
venção humana, fonte possível de comportamentos inespe- 
rados (ie. eventualmente não reconhecidos por outras 
componentes do sistema) é conseguida através do maior 
grau de flexibilidade incorporado. A contrapartida, como é 
habitual nos problemas de engenharia, é traduzida pela 
complexidade adicional exigida pela solução. Por questões 
de exposição, a apresentação cinge-se aos aspectos funda- 
mentais, deixando de lado pormenores de implementação. 
Por outro lado, o artigo foca apenas uma das subcompo- 
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(1) Características básica e específica são definidas no 

Cap. 3, ponto 17, de[1 ]. Aqui, é suficiente o exem- 

plo 1 apresentado na secção II.l, ponto 7. 
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nentes do projecto, sendo essencialmente extraído do capi- 
tulo 4 de [1]. Esta referência, e as referências nela incluí- 
das, detalham as outras subcomponentes do trabalho. 


I — INTRODUÇÃO 


Referir-se-á por controlo da produção, a actividade que 
acompanha a história da montagem dos produtos na linha 
de produção em série. 

Esta actividade tem o seu início quando é aberta uma 
obra de produção, que iremos denominar por lote. O con- 
trolo da produção é descrito de fotma sequencial no proce- 
dimento apresentado em pseudo-código na figura 1. 


procedimento CONTROLO DA PRODUCAO; 
begin 
identificação do lote; 
requisição de material; 
fornecimentos efectivos; 
while material em falta and lote aberto do 
begin 
requisição de material em falta; 
fornecimentos efectivos; 
end; 
end; 


FIGURA 1 — Controlo da produção 


Identificam-se os seguintes passos: 


1) A identificação do lote define os volumes de produção 
a serem montados de produtos com determinadas ca- 
racterísticas básica e específica (1). 


1) Com estes dados, procede-se à requisição de material. 
Com base na informação que define a composição dos 
produtos, é seleccionado o material a requisitar. Por 
análise das existências disponíveis, é feita a selecção das 
classes possíveis a serem requisitadas nos diferentes 
grupos de alternativas (ver 3.3, ponto 25 de [I]) (ver 
ainda secção IV mais abaixo). 

É feita uma previsão de material requisitado que irá ser 
fornecido. Todo o material requisitado é cativado para 
o lote em questão. 


lil) Os fornecimentos efectivos actualizam a informação de 
material previsto fornecido e permitem determinar se 
ainda há material requisitado não fornecido referido 
por material em falta. São criados movimentos de saída 
de material efectivamente fornecido. 


iv) Enquanto o material está em falta e o lote se mantém 


aberto, são desencadeadas sucessivas requisições de 
material em falta, seguidas da necessária informação 
dos fornecimentos que efectivamente abasteceram a li- 
nha de montagem. 


v) O processo termina quando o lote é completado ou 
quando é introduzida uma data de fecho para o lote. 


O procedimento da figura | descreve assim de um modo 
geral o encadeamento de actividades inerentes à história 
de um lote em particular. 

Estas actividades são resultado da comunicação entre 
três sistemas, nomeadamente os sistemas INFORMAT, 
(Informática), PRODUÇÃO e ARMAZÉM, 

O estudo do controlo da produção, que é apresentado 
neste artigo, está dividido em seis secções. As secções II, IH 
e IV referem-se ao processo individual de cada lote, en- 
quanto a secção V analisa as implicações da coexistência de 
vários lotes simultâneos no tempo. À secção II descreve a 
comunicação entre os sistemas e respectivas actividades. A 
secção Ill aprofunda as actividades locais e os processos de 
comunicação. A secção IV descreve o algoritmo de selecção 
de alternativas. Finalmente, a secção VI sumariza as con- 
clusões do artigo. 


II — COMUNICAÇÃO ENTRE SISTEMAS 
1. Introdução 


O procedimento simplificado da figura | apresenta as 
actividades, independentemente de qual o sistema que as 
executa. Procurou-se com esse resumo, definir a semântica 
do processo em si. Como sistemas intervenientes, e to- 
mando uma atitude simplificadora, assumir-se-ão apenas 
os sistemas INFORMAT, PRODUÇÃO e ARMAZÉM. A 
interacção entre estes sistemas no âmbito do controlo da 
produção será estudada nesta secção. 


2. Diagrama temporal de actividades e comunicação 
A cronologia de comunicação entre os três sistemas em 


causa, descrita informalmente na secção 1, é esquematizada 
na figura 2. 


São postas em evidência as actividades locais de cada 
sistema e a comunicação de informação entre sistemas. 


Notação: No diagrama de actividades e comunicação entre 
sistemas, usa-se a seguinte simbologia: 


Rs actividade local do sistema À 


CO >» — mensagem com origem no sistema A e 
destino o sistema B 
ES) N 
— — teste à informação local 
aaa — instrução interna ao esquema do 


diagrama 


+ 


INFORMAT PRODUÇÃO ARMAZEM 


requisição 


colectiva 


fornecimentos 
eieciivos 


lecho 
do lote 


FIGURA 2 — Diagrama temporal de actividades e 
comunicação 


3. Análise temporal do diagrama 
Cada eixo de tempo diz respeito a um sistema em parti- 


cular. 


Facto 1: Os eixos de tempo representados traduzem esca- 
las de tempo distintas para cada sistema. 


Por outro lado, o processo é iterativo. As iterações 
estão numeradas na variável i do diagrama e o bloco itera- 
tivo está assinalado devidamente. 


Facto 2: O processo exibe estrutura iterativa. O bloco ite- 
rativo tem os instantes de tempo significativos 
assinalados como função da respectiva iteração i. 


A precedência entre actividades e comunicação entre 


sistemas é expressa pela numeração adequada dos instan- 
tes de tempo significativos. 


Facto 3: Fora do bloco iterativo, verifica-se 
k> ko => W>tVkke No 
Dentro do bloco iterativo verifica-se 


DkKDk >>>). VkkeNVieN 
Di>Dji =>>wuD)>uD.VijeNVkeN 


lagli | 


Os instantes de tempo assinalados dizem respeito à comu- 
nicação entre sistemas. Assim, por exemplo, 


t; — Instante de recepção pelo sistema INFORMAT da 
mensagem de identificação de lote, emitida pelos sis- 
tema PRODUÇÃO. 


to — instante de envio ao sistema PRODUÇÃO do mapa 
de requisição colectiva, proveniente do sistema 
INFORMAT. 


Finalmente, o instante de tempo que assinala o fim da 
actividade é referenciado em cada um dos sistemas IN- 
FORMAT, PRODUÇÃO e ARMAZÉM por ta; (i), tos (i), 
€ tos (1), respectivamente. 


4. Análise de comunicação do diagrama 


O diagrama da figura 2 indica que o sistema ARMA- 
ZÉM não comunica directamente com o sistema INFOR- 
MAT no controlo da produção. O sistema PRODUÇÃO 
actua como meio de comunicação entre os sistemas IN- 
FORMAT e ARMAZÉM, Este facto coloca no sistema 
que recebe, e não no sistema que fornece, a responsabili- 
dade de controlar as entradas de dados (fornecimentos 
efectivos) no sistema INFORMAT. Porque é previsível, 
em futuro próximo, a inversão desta responsabilidade, com 
comunicação directa entre os sistemas ARMAZÉM e IN- 
FORMAT, este facto não afecta a implementação do con- 
trolo da produção. 


5. Análise de actividades do diagrama 


No decorrer do capítulo, são descritas em pormenor as 
actividades do sistema INFORMAT no que respeita ao 
controlo da produção. As actividades do sistema PRO- 
DUÇÃO são as de gestão do material que serve a monta- 
gem dos produtos de um lote. Do diagrama da figura 2 
concluem-se os seguintes pontos gerais, que caracterizam a 
função de cada sistema no processo de controlo da produ- 
ção. 


Facto 4: O sistema INFORMAT tem a seu cargo, o pro- 
cessamento automático de dados. 

Facto 5: O sistema PRODUÇÃO garante que há, no sis- 
tema ARMAZÉM, requisições de material, para 
que este seja abastecido à linha de montagem. 

Facto 6: O sistema ARMAZÉM abastece a linha de mon- 
tagem com o material que lhe é requisitado. 


II — ANÁLISE DAS ACTIVIDADES LOCAIS 
E MENSAGENS DE COMUNICAÇÃO 


KI.1. Identificação do lote 
6. Especificações 

Um lote é uma obra de produção. Em relação à identifi- 
cação do lote, enumeram-se as seguintes especificações, que 


estão subjacentes ao funcionamento da empresa: 


|) Um lote diz respeito a produtos com a mesma caracte- 
rística básica e diferentes características específicas. 


Para a montagem de produtos com características bási- 
cas distintas, são abertos lotes distintos. 

1) Um lote estabelece um volume de produção total dos 
produtos a montar. Este volume de produção é discri- 
minado por caracteristica específica. 


7. Consequências das especificações 


A especificação i) tem como consequência a existência 
em simultâneo no tempo, de vários lotes em montagem. 
Esta situação é analisada em detalhe na secção V. 

O exemplo seguinte ilustra o impacto da especifica- 


ção ii). 


Exemplo 1: Identifique-se a característica básica com um 
modelo de produtos e a característica específica com uma 
marca. Pela definição de características básica e específica 
(ver capítulo 3, ponto 17 de [1]), um modelo pode ter 
várias marcas. Seja o modelo MOD, e as m marcas res- 
pectivas MAR, MAR pj, MAR as. MAR ym- O lote 
L001 estabelece para o modelo MOD, o volume de produ- 
ção VPROD, = 1000. Se VPRODs,; representar, para este 
lote, o volume de produção a montar da marca MAR, 
verifica-se 


| VPROD,; = VPROD, 


im 3 


O lote fica totalmente identificado por 


— número do lote, (Ex: LOO1). 

— característica básica, (Ex: modelo MOD,). 

— volume de produção total, (Ex: VPROD, = 1000). 

— volume de produção discriminado por característica 
específica, (Ex: VPROD,;, i= 1,..., m). 


Esta informação é enviada pelo sistema PRODUÇÃO ao 
sistema INFORMAT, como se pode ver na figura 2. A 
actividade de requisição de material é desencadeada no sis- 
tema INFORMA, após validação da identificação do lote 
recebida. 


H1.2. Requisição de material 
8. Entrada e saída de dados 


A requisição de material é uma actividade do sistema 
INFORMAT, que tem como entrada de dados: 


— a estrutura de dados de incorporação de items nos 
produtos (relações INCORB e INCORE, global- 
mente denominadas INCOR), 

— os dados de identificação do lote, 

— o nível de existências dos items (relação EXIST), 

— a descrição dos items (relação DESCRIT). 


Como saida de dados é criada a relação de quantidades 
de items incorporados a requisitar para a montagem do lote 
(relação REQLOTE); ver o esquema de blocos na figura 3. 
Por razões de espaço, não é aqui possível detalhar o mo- 
delo lógico de dados, nem a estrutura física de suporte. 
Para fácil compreensão do artigo, as relações anteriores 
podem ser assimiladas pelo leitor de forma intuitiva a tabe- 


las ou ficheiros convencionais. Assim, por exemplo, a rela- 
ção EXIST definirá para cada item qual a sua existência em 
armazém. 


! Requisição de 
Material 


REQLOTE (CODIT, DESIGIT, 
ARM, CCUST, QREQ) 


FIGURA 3 — Diagrama de entrada saida da actividade 
de requisição de material, 


Na figura 4 apresenta-se o procedimento de requisição de 
material, 


procedimento REQUISICAO DE MATERIAL; 
begin 
requisição de itens nao pertencentes 
aos grupos de alternativas CREQLOTE |); 
selecção e clusses possíveis a requisitar 
para cada grupo de arternativas CREQLOTE , ); 
REQLOTE = REQLOTE, U REQLOTE 
end; 


É 


FIGURA 4 — Requisição de material 


Para o item 1,, que não é elemento de nenhum grupo de 
alternativas, i.e. 1; e Go, a quantidade a requisitar Qreq é 
função da quantidade incorporada Q,.. e do volume de 
produção do lote. 

Seja o lote k respeitante a produtos com característica 
básica CBAS,. O volume de produção VPROD, do lote k 
está discriminado por característica específica CESP;. Veri- 
fica-se, como já se viu em II.l, 


m 
£ VPROD,; = VPROD,, 


i=1 


em que m é o número total de características específicas 
(p.e. marcas). A quantidade Q,eg é obtida por 


m 
Ve Go, QreglD= E Qinc (1, CBAS,, CESP;) * VPROD;; 


a 
ss: 
| — 


Exemplo 2: Seja o lote LO0I com a identificação introdu- 
zida no exemplo 1. Seja REQLOTE ,,, a relação resultante 
do processo de requisição de items não pertencentes aos 
grupos de alternativas. Na figura 5, REQLOTE, é descrita 
em linguagem pseudo-código. Os símbolos m, q e |X| da 
figura 5 têm o significado habitual dos operadores projec- 
ção, selecção e junção natural da álgebra relacional [2]. 


INCORLOT = TT (oe (INCOR)) 
CESP, CODIT, CCUST CHAS=MOD, ão 
PARTE ,QUANT, UNID a GRUP=0 4 QUANT fu 
a REQ =1 
while not eof INCORLOT 
begin 
read INCORLOT (CESP, CODIT, CCUST, FANIE, QUANT, UNID); 
for all à, 
if CESP=MAR | then QREQ:= QUANT * VPROD,/UN4L 
write LHCORLOT (CESP, CODIT, CCUST, PARTE,QREO, UNID); 


end; 
REQLOTE = = (m (IHCORLOT) |x| DESCRIT) 
CODIT, CODIT, 
DESIGIT, CCUST, 
ARM, QREy 
CCUST, 
QREO 


FIGURA 5 — Definição em pseudo-código da relação REQLOTE ,. 


REQLOTE, é a relação de requisição de items pertencentes 
a grupos de alternativas. Verifica-se REQLOTE = REQ- 
LOTE, U REQLOTE,. 

É adiada para a secção IV a apresentação detalhada do 
algoritmo de selecção de alternativas. 


9, Previsão de fornecimentos 


Para a emissão dos movimentos de saida de material do 
ARMAZÉM, o sistema INFORMAT aguarda até ao ins- 
tante tg (ver figura 2), em que lhe é dada informação sobre 
os fornecimentos efectivos de cada item. 

Introduza-se, para cada item, o atributo de existência 
disponível EXSDISP, distinto da existência real do item, 
EXS. No ponto 10, é definido exactamente o valor de 
EXSDISP para cada item. De momento, admita-se que a 
existência disponível é uma fracção da existência real de 
cada item. Veja-se a tabela | para a especificação da no- 
menclatura usada no seguimento da apresentação. 


EXS - EXiStência real 

EXSDISP = Existência  DISPonivel 

QREQ - Quantidade REQuisitada 

QFORP  -  Quantigaue ruxnecida Prevista 
QFOR - Quantidade FURnecida efectivamente 
QCAT - Quantidade CATiva 


TABELA 1 — Variáveis de processamento da requisição 
de material, 


Facto 7: A previsão de fornecimentos é obtida do seguinte 
modo: 


if EXSDISP < QREQ 


then QFORP = EXSDISP 
else QFORP = QREQ. 


Estuda-se na subsecção 111.4 que os fornecimentos efec- 
tivos apenas terão que ser comunicados, caso contrariem 


esta regra. Isto é, para reduzir o fluxo de comunicação, 
adoptar-se-á o tratamento explícito apenas das excepções. 


10, Existência disponível — Relação EXCAT. 


Para cada item, a quantidade requisitada traduz a quan- 
tidade efectivamente necessária para a montagem dos pro- 
dutos em causa, já que se baseia nos dados da relação 
INCOR. A necessidade de prevenir a requisição de material 
ao armazém no lapso de tempo ty — ts > O (ver figura 2) 
sem com isso perturbar o controlo da produção leva à in- 
trodução do conceito de quantidade cativa. 

A requisição de material cativa para o lote, cujo processo 
está a decorrer, as quantidades requisitadas para cada item. 

Com o controlo da produção, surge assim uma outra 
relação básica para o sistema INFORMAT, onde se man- 
tém a informação das existências cativas para cada item, 

Seja a relação EXCAT (EXistências CA Tivas) com atri- 
butos CODIT, CCUST, OBRA, QCAT e dependências 
funcionais DFEC. 


DFEC = | CODIT: CCUST: OBRA  -s QCAT) 


Resultado 1: CODIT: CCUST: OBRA é a única chave da 
relação EXCAT, e esta encontra-se na forma 
normalizada de Boyce-Codd. 


Definição 1: A existência disponivel EXSDISP, para o 
item A, obtem-se a partir da sua existência real, 
EXS e quantidade cativas, pelo seguinte pro- 
cedimento: 


= 


EXSDISP=EXS 
while not eof EXCAT 
begin 
read EXCAT (CODIT, CCUST, Ubna,QCAT) 
if CODIT=A 
then EXSDISP=EXSDISP=QCAT; 
end; 


Está subjacente na relação EXCAT a existência de várias 
obras LOTE simultâneas. Este ponto vai ser discutido na 
secção V. 


1.3. Mapa de requisição colectiva 


O mapa de requisição colectiva é a mensagem enviada 
pelo sistema INFORMAT ao sistema PRODUÇÃO. É 
uma enumeração exaustiva dos items e respectivas quanti- 
dades necessárias para a montagem de determinados pro- 
dutos. 


11. Atributos da mensagem 


Uma mensagem é apresentada neste estudo, como uma 
relação (ver capítulo 2 de [1]). O mapa de requisição colec- 
tiva tem os campos ou atributos apresentados na tabela 2. 


CODIT  -  CODigo do Llem 

DESIGIT - DESIGnação do ITem 

ARH - ARMaz &m 

EXSDIN -  EXiStência Lisponivel INicial 
(antes da presente requisição) 

QJAFOR —-  Quantiaaoe JA FORnecida 

QREQ - Quantidade REQuisitada 

QFOR - Quantidade Fuknevida efectivamente 


TABELA 2 — Atributos do mapa de requisição colectiva 


Em [3] são apresentados exemplos concretos de mapas 
de requisição colectiva. 


12. Cabeçalho do mapa de requisição colectiva 


O mapa de requisição colectiva é um documento de mo- 
vimentação de items, em que as quantidades movimentadas 
são as quantidades fornecidas efectivamente. 

O cabeçalho do mapa tem os atributos numerados na 
tabela 3, comuns a todos os tuplos do mapa, 


NUMD - NiMero do Documento 

OBRA - OBRA (lote) 

CBAS - Caracteristicas  BASica 
VPROD - Volume de PHODução total 


par Característica ESPecifica, Volume 


i 


(CESP, VPROD, ) 
: de PRODução para i=l,...,N, em que m 
ê o número de características a que à 


obra diz respeito 


TABELA 3 — Atributos do cabeçalho do mapa de 
requisição colectiva. 


Como se conclui do diagrama da figura 2, o mapa de 
requisição colectiva é uma mensagem de comunicação en- 
tre os três sistemas INFORMAT, PRODUÇÃO e AR- 
MAZÉM. Estuda-se a acção de cada sistema em relação a 
este mapa. 


13. Acção do sistema INFORMAT 


O mapa de requisição colectiva é gerado automatica- 
mente pelo sistema INFORMAT com base em: 


1) relação REQLOTE de saída do processo de requisição 
de material (ver figura 3); 
1) informação residente sobre a existência disponível. 


A partir dos dados de identificação do lote, o sistema 
INFORMAT preenche os atributos OBRA, CBAS, VPROD, 
CESP, e VPROD, (i= 1,..., m) do cabeçalho. Quanto ao 
atributo NUMD, o seu valor é gerado automaticamente 
por incremento do número do último documento emitido. 

Os atributos do mapa, ver tabela 2, e à excepção de 


QFOR, vêm directamente da relação REQLOTE. A exis- 
tência disponível inicial EXSDIN, é obtida pelo processa- 
mento que define a existência disponível (ver definição 1), 
no instante inicial da requisição. 

O atributo QFOR não é preenchido pelo sistema IN- 
FORMAT. Será o sistema ARMAZÉM que o preencherá 
informando das quantidades efectivamente fornecidas. 

A situação em que EXSDIN < QREQ é assinalada com 
um alerta na existência disponível inicial (ver [3]). 


14. Acção do sistema PRODUÇÃO 


Como se vê na figura 2, o sistema PRODUÇÃO faz uma 
correcção à requisição. 


Facto 8: O sistema PRODUÇÃO tem autoridade para al- 
terar o conteúdo do campo da quantidade requisi- 
tada, QREQ. 


Relembre-se que (ver figura 5) 
QREQ = VPROD * INCOR. QUANT/INCOR. UNID 


A alteração da quantidade requisitada QREQ é conse- 
quência das seguintes situações possíveis, detectadas pelo 
sistema PRODUÇÃO: 


|) Os valores de INCOR. QUANT ou INCOR, UNID não 
estão actualizados. 


1) Por informação da linha de montagem, o sistema PRO- 
DUÇÃO tem conhecimento de que há um excesso ou uma 
falta de material que justifica a diminuição ou o aumen- 
to da quantidade requisitada (por exemplo, resultante de 
peças defeituosas ou sobras de lotes anteriores acumula- 
das na linha). 


Enquanto a situação ii) não exige nenhuma acção ex- 
traordinária por parte do sistema PRODUÇÃO, para além 
da correcção justificada da quantidade requisitada pelo sis- 
tema INFORMAT, o mesmo não sucede em relação à si- 
tuação i). A situação 1) exige que o sistema PRODUÇÃO 
desencadeie uma acção de envio de mensagem ao sistema 
que mantém actualizada a relação INCOR, que é o sistema 
TÉCNICA (ver tabela 3.8 de [1]). informando-o da ano- 
malia detectada. 


Facto 9: O sistema PRODUÇÃO tem autoridade para 
acrescentar novas requisições de material, através 
da inserção dos atributos CODIT, DESIGIT, 
ARM, QREQ. 


Este facto traduz uma situação de que o facto 8 é um 
caso particular. 


15. Acção do sistema ARMAZÉM 


O sistema ARMAZÉM recebe o mapa de requisição co- 
lectiva e procede ao abastecimento dos items listados nas 


quantidades requisitadas corrigidas ou não pelo sistema 
PRODUÇÃO. 


Facto 10 O sistema ARMAZÉM introduz a informação 
quanto à quantidade efectivamente fornecida 
QFOR. 


Notação: A notação (x) refere-se ao valor do atributo x de 
uma mensagem. 


Facto 11 O ARMAZÉM atribui valores a QFOR, de 
acordo com a seguinte regra: 


Se a quantidade fornecida no abastecimento é igual a 


(QREQ) 
então (QFOR) = «/ sinal de acordo ) 
senão (QFOR) = valor da quantidade fornecida. 


Denomina-se por QFOR indistintamente a designação 
do atributo ou o seu conteúdo. O contexto onde é inserida 
a variável torna claro o sentido pretendido. 

Com respeito ao intervalo de tempo de abastecimento 
(acção do sistema ARMAZÉM), conclui-se da figura 2, o 
facto 12. 


Facto 12: t6 — 15 2> mel) — tsQ), re No 


Este facto traduz a constatação de que um abastecimento 
não é uma actividade com duração própria, antes depen- 
dendo do volume de items a abastecer. Como será estudado 
posteriormente, o volume de items requisitados no mapa de 
requisição colectiva é largamente superior ao volume de 
items requisitado numa requisição de faltas. 

Serão postas presentemente em evidência, duas situações 
distintas que conduzem ao não abastecimento do material 
requisitado na sua totalidade. 


16. Não abastecimento por insuficiência de 
existência física 


O sistema ARMAZÉM apenas pode fornecer quantida- 
des inferiores ou iguais às que existem fisicamente. 


17. Não abastecimento por outros condicionamentos 


Apesar da existência poder permitir o fornecimento inte- 
gral, há condicionamentos particulares que o impedem. Es- 
tão neste caso, os items cujo volume físico ou fragilidade 
não permita a sua acumulação na linha de montagem 
durante o periodo total de fabrico do lote. Nesta situação, é 
fornecida apenas uma fracção da quantidade requisitada, 

Todo o material não fornecido, fica em falta para o lote, 
sendo requisitado de novo em iterações seguintes. 


H1I.4. Fornecimentos efectivos 


18. Fornecimentos efectivos do mapa de requisição 
colectiva 


Depois de receber o mapa de requisição colectiva prove- 
niente do sistema ARMAZÉM (ver figura 2. ty > ty),o 
sistema PRODUÇÃO comunica ao sistema INFORMAT 
os fornecimentos efectivos abastecidos pelo ARMAZÉM. 
Estes são devidamente assinalados no atributo QFOR, de 
acordo com a regra enunciada no facto 11. 

O suporte para a comunicação tem os mesmos atributos 
que o mapa de requisição colectiva. Atendendo às quanti- 
dades previstas fornecidas (ver [1.2, ponto 9), só há neces- 
sidade de comunicar ao sistema INFORMAT a informação 
que as contraria — regra de excepção referida anteriormente 
e formalizada no facto 13. 


Notação: A notação (x)a refere-se ao valor do atributo x 
de uma mensagem, atribuído pelo sistema A. 


Facto 13: Osistema PRODUÇÃO apenas informa o sistema 
INFORMAT, nas seguintes situações: 


1. (QREQ)pronução * (QREQ)inrorMmaT 
2. (EXSDIN) 2 (QREQ) e (QFOR) * 
3. (EXSDIN) < (QREQ) e (QFOR)* (EXSDIN) 


A situação | traduz os factos 8 e 9. As situações 2 € 3 
resultam de má previsão, já que contrariam as regras apre- 
sentadas no facto 7. 

Na situação 2, a existência disponível inicial EXSDIN é 
superior ou igual à quantidade requisitada QREQ. A quan- 
tidade prevista fornecida QFORP é igual à quantidade re- 
quisitada, de acordo com o facto 7. O valor da quantidade 
efectivamente fornecida (QFOR) diferente de «/ significa 
que a quantidade requisitada QREQ não foi fornecida na 
totalidade (ver facto 11), o que contraria a previsão QFORP. 
Na situação 3, a quantidade prevista fornecida QFORP é 
igual à existência disponível inicial EXSDIN, pois esta é 
inferior à quantidade requisitada QREQ. A previsão é con- 
trariada pelo não fornecimento efectivo da existência dis- 
ponível EXSDIN. expresso em (QFOR) * (EXSDIN). 

Uma análise estatística da frequência de ocorrência, de 
cada uma destas situações, conduz à justificação da vali- 
dade da previsão de material fornecido. 


19. Análise estatística 


Para o tuplo genérico de um mapa genérico, definam-se 
os seguintes acontecimentos: 


A = ocorrência de uma das situações descritas no 
facto 13, 


B; = ocorrência da situação i particular. 


Porque os acontecimentos B; são mutuamente exclusivos, 


3 
Facto 14: Pr[A]) -=- E 


Seja N o número médio de tuplos de um mapa, L.e., N é 
o número médio de requisições individuais de items numa 
requisição colectiva de material. 


Facto 15: O número médio de tuplos N.,. exigindo a acção 
de excepção é 


Nexo = PrfAjJ.N 


Por análise estatística dos acontecimentos, resultaram os 
seguintes valores numéricos: 


0.010, i 
Pr (Bj = 40.040, 1= 
0.050, 1 


II 


I 
3 + 


donde Pr [ A] = 0.1. Este valor quantifica a economia 
de esforço representada pela regra de excepção adoptada, 


Pelo ponto 16, este valor será tanto menor, quanto mais 
fiável for a informação residente no sistema INFORMAT. 

Os valores de Pr [ B;] e Pr [| B;) justificam-se pelo 
decorrer simultâneo de vários lotes, como será focado na 
secção V. De facto, embora a requisição de material fun- 
cione sequencialmente, o abastecimento no sistema AR- 
MAZÉM não segue nenhuma ordem específica. Assim as 
existências disponíveis no acto de fornecimento podem não 
coincidir com as disponíveis EXSDIN, no momento da 
requisição de material. 

Quanto à Pr (B,], pode ainda decompôr-se o seu va- 
lor em duas parcelas de acordo com as alternativas 1) e ii) 
indicadas no ponto 14. Uma análise estatística mais fina 
mostrou que é a situação ii) — excesso ou falta de material 
na linha de montagem — que contribui mais significativa- 
mente para o valor de Pr [ B,). 


20. Fornecimentos efectivos da lista de requisição 
de faltas 


Por pedido de requisição de faltas com origem no sis- 
tema PRODUÇÃO e destino o sistema INFORMAT é ge- 
rada a lista de material ainda em falta para a montagem do 
lote. O sistema PRODUÇÃO envia esta lista ao sistema 
ARMAZÉM, que a preenche e devolve de novo ao sistema 
PRODUÇÃO. 

Sempre que é recebida a lista de requisição de faltas, o 
sistema PRODUÇÃO comunica ao sistema INFORMAT 
os fornecimentos efectivos. A regra seguida é distinta da 
que foi indicada no facto 13 para o mapa de requisição 
colectiva. A análise revela que os acontecimentos que ocor- 
rem são distintos dos discriminados anteriormente. Por 
outro lado, as respectivas possibilidades de ocorrência são 
também significativamente diferentes, não se justificando a 
regra de excepção. Esta constatação não é de estranhar, 
dado que se trata agora de acompanhar os fornecimentos 
das excepções detectadas aquando do mapa de requisição 
colectiva. 


Facto 16: O sistema PRODUÇÃO apenas informa o sistema 
INFORMA, nas seguintes situações: 


|. (QREQ)produção * (QREQ)inrorMaT 
2. (QFOR) * O 


Os valores elevados da probabilidade de ocorrência des- 
tas situações são ainda justificados pelo significado do 
mapa de requisição colectiva versus a lista de material em 
falta, e pelos respectivos tempos de abastecimento. 

Como foi posto em evidência no facto 12, e atendendo à 
figura 2, verifica-se que tg — ts >> tus (1) — tis (1), 
Le N. 


Facto 17: De modo a ser conhecido o abastecimento do 
material em falta com a assiduidade desejada, o 
intervalo de tempo ty (1) — tis (), ie N, é 
muito menor que tg — ts. 


Este facto, apoiado numa análise estatística que estuda a 
frequência das situações descritas no facto 16, justifica a 
regra estabelecida de necessidade de comunicação. 

Para concluir o estudo da comunicação de fornecimentos 
efectivos, faz-se notar que a regra de excepção adoptada 
procurou satisfazer um princípio geral de acção, que subjaz 
a filosofia de processamentos distribuídos. 


IO 


Princípio: Em processamentos envolvendo sistemas distin- 
tos é desejável minimizar o volume de dados de 
comunicação entre sistemas. 


A aplicação deste princípio, ao caso particular em es- 
tudo, levou não só ao processo de previsão de material 
fornecido, como ainda à análise estatística e semântica das 
diferentes situações que poderiam conduzir à necessária 
comunicação de dados. 


HI.5. Actualização de dados e movimentos de saída 


A actualização de informação local no sistema INFOR- 
MAT é completada com a comunicação dos fornecimentos 
efectivos, que representam saída de material do sistema 
ARMAZÉM. 

No capítulo 3, subsecção 3.2.5, ponto 12 de [1], é apre- 
sentada a relação MOVIT que traduz a movimentação dos 
items, que é actualizada pela actividade em estudo do con- 
trolo da produção. Por outro lado, é automaticamente ca- 
tivado o material que se mantém em falta. 

O acesso simultâneo por parte de diferentes sistemas à 
relação MOVIT é estudado no capítulo 7 de [1], onde é 
desenvolvido um algoritmo para controlo de concorrência 
de acessos. Razões de espaço levam à omissão destes tópi- 
cos no presente artigo. 


1.6. Requisição de material em falta e fecho do lote 


O sistema PRODUÇÃO emite um pedido de requisição 
de faltas, identificando o respectivo lote. Depois de testar se 
o lote em causa está aberto, o sistema INFORMAT pro- 
cede à actividade de requisição do material em falta. 


Facto 18: É requisitado o material que se encontra ainda 
cativo para este lote. 


O processo de requisição de faltas é iterativo. 
A actualização de dados mantém a iteração k depen- 
dente dos dados das iterações anteriores. 


Notação: X* representa o valor do atributo X na iteração k. 


Seja k a presente iteração. Então, tendo em conta as 
abreviaturas apresentadas na tabela 2, secção 111.3, ponto LI, 


QREQ* = QREQE—! — QFOR*k-! 


QIAFOR* - >,  QFOR' 


Já se justificou a existência de material em falta (ver pontos 
l6 e 17, secção 11.3). 

A quantidade já fornecida QJAFOR é um acumulado 
das quantidades efectivamente fornecidas nas iterações an- 
teriores. A relação entre QREQ* e QJAFORÉ permite ter 
noção da fracção de material ainda em falta. 

A acção de fecho de lote pode ter origem em duas situa- 
ções não exclusivas: 


1) abastecimento integral do material em falta, 
ii) decisão do sistema PRODUÇÃO. 


Na situação 1), se o material não foi integralmente 
abastecido, há uma actualização do material cativo, liber- 
tando o material em falta do seu estado de cativo para o 
lote fechado. 


IV. ALGORITMO DE SELECÇÃO DE ALTERNATIVAS 
21. Introdução 


No capítulo 3 de [1] é apresentada a constituição dum 
produto como sujeita à definição de alternativas. Para a 
requisição de material com vista à montagem de produtos 
específicos, torna-se necessário estabelecer um algoritmo de 
selecção de alternativas. Subconjuntos de items (não neces- 
sariamente disjuntos) (com a respectiva incorporação defi- 
nida) são ditos em alternativa se a constituição de um pro- 
duto implicar a integração de um e um só desses 
subconjuntos. 

No referido capítulo 3, distinguem-se dois casos depen- 
dendo do critério adoptado para a selecção de material por 
grupo de alternativas; 


CRITÉRIO 1: selecção por ordem de prioridades 
pré-estabelecida. 


CRITÉRIO 2: selecção por busca de stock mínimo. 


O caso | será abreviadamente referido como critério de 
prioridades. enquanto que o caso 2 será referido como o 
critério da busca do stock minimo. O stock minimo não é a 
existência mínima, mas sim esta normalizada pela incorpo- 
ração. 

Um grupo de alternativas verifica o critério 2, se não foi 
estabelecida uma ordem de prioridades de escolha entre as 
classes possíveis que o compõem. A busca do stock minimo 
consiste na procura do item do grupo de alternativas, cujo 
stock seja o minimo. O item de stock minimo determina a 
primeira escolha de uma classe possível do grupo a que 
pertence. Se o stock mínimo não permitir atingir o volume 
da produção total do lote, nova busca de stock mínimo de 
entre os items das restantes classes possíveis, determina a 
segunda classe escolhida. A busca iterativa termina, 
quando se conseguiu atingir o volume de produção total do 
lote ou se esgotarem as classes possíveis. 

Nos dois casos. se se esgotarem as classes possíveis, sem 
ser atingido o volume de produção total, o algoritmo requi- 
sita o volume de produção em falta, de uma classe possível 
aleatória ou pré-estabelecida. 

No decorrer desta secção, estuda-se em detalhe cada um 
dos casos possíveis, bem como o procedimento algoritmico 
implementado. A apresentação é acompanhada de exem- 
plos ilustrativos. 


IV.l. Selecção por ordem de prioridades pré-estabelecida 


O algoritmo tem como entrada os dados sobre a incorpo- 
ração dos items nas diferentes classes possíveis e determina 
os volumes de requisição para os items de cada classe pos- 
sivel. 


Entrada: |. G é um grupo de alternativas. 
2. nel é o número de classes possíveis para o gru- 
po G. 
3. PRIOR [ncl] é a matriz de prioridades por 
classe possível. PRIOR [i] =) significa que a 
classe | tem prioridade 1. 


| 


4. nit é o número de items do grupo de alternati- 
vas. 

5. QINC [nit x nel] é a matriz de quantidades 
de incorporação por item e por classe possível. 

6. EXSD [nit] é o vector de existências dispo- 
níveis para cada item, 

7. VPROD é o volume de produção a atingir. 


Objectivo: Determinar a matriz VCLAS [ncl] dos volu- 
mes a requisitar de cada classe possível de modo 
ao btera matriz QITEM [nit] de quantidades 
a requisitar de cada item do grupo G, tal que 


|. QINC x VCLAS = QITEM 


nel 
2. z VCLAS [k] = VPROD. 
k = 1 
3. Se obedeceu às prioridades definidas em 
PRIOR. 


O algoritmo é explicado com o auxílio de um exemplo 
particular. O procedimento que implementa este algoritmo 
é apresentado na figura 7. O código escolhido foi o Pascal, 
Já que se considera esta uma linguagem auto explicativa. 


22. Exemplo 3 


Seja o grupo de alternativas representado na figura 6 
com 7 items associados em 3 classes possíveis. 


FIGURA 6 — Items e respectivas classes possíveis 


Admitam-se a matriz de incorporações QINC e a matriz 
de existências disponíveis EXSD: 


QINC = 


EXSD = [450 825 1000 125 800 632 2500]! 


O volume de produção assumido é VPROD = 1000 con- 
siderando-se o vector de prioridades PRIOR = [23 1]T. 
As classes 1, 2 e 3 têm respectivamente prioridades 3, 1 e 2. 

Na iteração | do algoritmo é seleccionada a classe de 
prioridade mais elevada, que no exemplo em consideração 


é a classe 2. Começando na classe 2, determina-se o volume 
máximo permitido pelas existências disponíveis. Os ele- 
mentos da classe 2 são 14, 1, e 1y, com existências dispo- 
níveis, 825, 1000 e 2500. 

O volume a requisitar da classe 2, será: 


825 1000 2500 
VCLAS [2] = L min ( | Tr VPROD )| 


onde L call representa a parte inteira do argumento e min 
(aj. à; ... an) calcula o minimo entre aj, à2 ..., ap. Os 
argumentos da função min correspondem aos valores das 
existências disponíveis divididos pela quantidade de incor- 
poração dos items da classe 2. Para o exemplo em estudo, 
com VPROD = 1000, 


VCLAS [2] = 500 
como 
VCLAS [2] < VPROD = 1000, 


a classe de prioridade seguinte é analisada, depois de o 
algoritmo actualizar o vector das existências disponíveis 


EXSD = [455 325 0 125 800 632 1000]T, 
e o volume de produção, 
VPROD = 1000 — VCLAS [2] = 500. 


Para a classe de prioridade 2, ou seja a classe 3, obtém-se 
de forma idêntica, 


vciaspj=lLmint 4 — — — vpRrOD)I 
3 | 


O item 1, é comum às classes 2 e 3. Porque foi este o item 
que já decidiu o volume de produção possivel na iteração 
anterior, isto é da classe 2, a sua existência disponível foi 
completamente esgotada EXSD [3] = 0. O volume de 
produção possível com a classe 3 é consequentemente nulo. 
O algoritmo passa imediatamente à classe seguinte, neste 
caso à última classe possível, ou seja, a classe |, de priori- 
dade 3. 
O volume de produção possivel vem 


455 325 BOO 632 1000 


2 | 2 3 | 


vcLAS [1]= Lmin ( VPROD) À. 


Como VPROD = 500, vem 
VCLAS [1] = 210. 


As existências disponíveis são uma vez mais actualiza- 
das. O volume de produção total não foi atingido 


VCLAS [2] + VCLAS [1] =710 < 1000, 
Para suprir o débito 1000 — 710 = 290, admita-se que a 
classe possivel a ser requisitada nesta situação é pré-estabe- 


lecida como sendo a de prioridade máxima, neste caso, a 
classe 2. É actualizado assim o valor de VCLAS [21 


VCLAS [2] = 500 + 290 = 790. 


1 


12 


As existências disponíveis são novamente actualizadas de 
acordo com esta nova requisição. Note-se no entanto que 
as existências disponíveis são limitadas inferiormente ao 
valor zero. O algoritmo, neste exemplo particular, teria 
como saida as seguintes matrizes VCLAS de volumes de 
produção por classe e QITEM de quantidades requisitadas 
por classe. 


VCLAS =[210 790 0]! 
QITEM = [420 1000 1580 O 420 630 2580]" 
A existência disponível para cada item é dada por 


EXSD = [30 00 125 380 2 0]!. 


procedure SEL PRIOR; [sulceção por prioridades! 


cunst 
NiEs... Intmero de items do grupo) 
nelas... intmero de classes pusslveis] 
infinitos... 

var 


À, ), P. pmax, pcl: integer; 
VIAPOS, VPROU: reui; 

FLM: boolean; 

EXSD, QITEM: array [1..nit] of real; 
VCLAS: array [1..nci] of real; 
PRIOR: array [1..mel] of integer; 
QINC: array [1..nit,l..ncl] of real; 


function min (x, y : real): real; 
(determina o minimu entre x ey] 


begin 
[inicialização de variáveis] 
for 1:=1 to nel do, 
regu (Fniun [1]); 
for 1:=1 to nit do 
read (EXSD [1]); 
for 1:=1 to nit do 
fur j:=1 to nel do 
reao (QINC [1,7]; 
for 1:=1 tó nit do 
QITEM [1):=0; 
for 1:=! To nel do 
VCLAS [il:=infinito: 
VJ APOS: =O 
FIM: =false; 
pmax:=1; [priorioaoe máxima & 1] 


[cálculo de VCAS) 
while bt (FIM or p=ncl) do 
DegIN 
prspel; 
pel:=PRIOR Lp); [a classe pcl tem prioridade p] 
for 1:=1 to rut 00 
1f not (QINC [i,pcl]=0) 
then VCLAS [pclJ:="min (VCLAS [pel], 
[(Lasu La J/QINC [4,pcl])]; 
if not (VCLAS [pcl]=0) 


then begin 
if VOLAS [pelj > vrnUv-U APOS 
then begin 
VELAS [pcl]:=VPROD-VJ APOS; 
FIM: =true; 


end; tnen, df] 
WIAPOS:=VJ APOS+VCLAS [pel]; 
for 1:=1 to nit do 
EXSD [i):=Ex5u |1)=ViLAS [pcl) * 
QINCIi,pcll; 


end; then, df) 
end; while) 
for i:=) ly mui do 
if VCLAS [ilsinfinito then VCLAS [1]:=0; 
if not (VPROD=V APOS) 
Laser begin 
VCLAS [PRIOR [pmax]]:=VCLAS [PRIOR [pmax]) 
+ (VPROD=-VI Abusos) 
For 1:=1 to nit do 
àf not (QINC [i, PRIOR [pmax)]=0 
then begin 
EXSD L1J:=tASU [(4]-(VEROL-V] APOS) 4 
E QINC [i, PRIOR [quiax]]; 
if EXSD [1] < O 
tneri EASU [1]:=0: 
Ena; 


end; 


Tor 1;=1 to nit do [requísição de materius] 


fer j:=1 to mel do 
QITEM [1]:=QITEM [4] + VELAS [4] * QUIC [1,3]; 
er, 


FIGURA 7— Algoritmo de selecção de alternativas com 
prioridades pré-estabelecidas 


O procedimento apresentado na figura 7 é uma versão sim- 
plificada do algoritmo implementado. Pretende-se apenas 
realçar o aspecto da selecção das classes possíveis, não se 
evidenciando pormenores de acesso e actualização de da- 
dos nas relações respectivas. Daí que a inicialização de da- 
dos consista na entrada de matrizes. Na implementação 
prática, esta entrada é equivalente à leitura dos dados nas 
respectivas relações. As instruções de read devem ser assim 
interpretadas, como procedimentos de leitura dos dados 
nas respectivas relações, com criação da estrutura de dados 
utilizada pelo algoritmo. O tempo de processamento da 
busca iterativa é minimizado pela manipulação de dados 
estruturados matricialmente, pois decorre em memória cen- 
tral, São assim minimizados os acessos à memória de 
massa, para optimização do tempo de resposta, 


IV.2. Selecção por busca do stock mínimo 


O algoritmo, neste caso, procede de modo distinto, embora 
tenha como base os mesmos dados. 


Entrada: 1. G é um grupo de alternativas. 

2. nel é o número de classe possíveis para o gru- 
po 6. 

3. nit é o número de items do grupo de alterna- 
tivas. 

4. QINC [nit x nclJé a matriz de quantidades 
de incorporação por item e por classe possível. 

5. EXSD [nit Jé o vector de existências dispo- 
níveis para cada item. 

6. VPROD é o volume de produção total. 


Objectivo: Determinar a matriz VCLAS [ncl ] dos volu- 
mes a requisitar de cada classe possível de modo 
a obter a matriz QITEM [nit ] de quantidades 
a requisitar de cada item do grupo G, tal que 


|. QINC x VCLAS = QITEM 


nel 
2. É, VCLAS [k] = VPROD. 
3. Se obedeceu à busca de stock mínimo. 


O algoritmo é explicado por recurso ao mesmo exemplo 
da subsecção IV.l. Na figura 8 é apresentado um procedi- 
mento em Pascal para a implementação do algoritmo. 

Considere-se o grupo de alternativas mtroduzido na fi- 
gura 6, com matriz de quantidades incorporadas QINC e 
existências disponíveis EXSD, dadas no ponto 22, da 
mesma subsecção IV.l. 

Seja mais uma vez VPROD = 1000. O algoritmo baseia- 
-se na construção da matriz EXSREST [nit x ncl]. Esta 
matriz traduz as existências para cada item que estariam 
disponíveis, se para o VPROD em causa, o item fosse re- 
quisitado para dada classe possível. 

Se para dada classe, o item tem incorporação nula, a 
respectiva entrada é inicializada a infinito. Para o grupo em 
exemplo, a matriz EXSREST vem 


EXSREST = | —I550 — 175 a o —I200 —2368 1500 
o — 175 — 000 “o so o —S0 | 
o — 1175 —2000 — 2875 —200 —368 ” 


em que, por exemplo, 


EXREST (6,1) = EXSD (6) — QINC (6,1) * VPROD 
= 632 — 3 * 1000 
= 2368 


Por busca do elemento mínimo da matriz EXSREST, o 
algoritmo identifica o item que conduz ao stock mínimo. A 
classe possível a que esse item pertence é a primeira a ser 
requisitada pelo algoritmo. Em caso do item de stock mi- 
nimo pertencer a mais que uma classe, é arbitrariamente 
escolhida a classe de menor número. O processo é iterativo 
e termina quando se atingiu o volume de produção dese- 
jado, ou quando já não há mais classes possíveis a requisi- 
tar. 

No exemplo em estudo, o elemento minimo da matriz é 
EXSREST [4,3 ] = — 2875. Na primeira iteração, é requi- 
sitada a classe 3. O volume de produção máximo permitido 
para esta classe possivel é determinado por processo idêntico 
à selecção por prioridades. 

Assim, os elementos da classe possível 3são 15, 14, 14, Is 
e Ig com existências disponíveis respectivamente 825, 1000, 
125. 800, 632. O volume de produção permitido para a 
classe 3 é 


| 825 1000 125 800 632 
vCLAs)= Lmin( — — =, —, —, vrrOD) À, 
2 3 3 Tt 


Com VPROD = 1000, vem VCLAS [3 ]=41. 


O algoritmo actualiza o vector das existências disponíveis. 
O volume de produção em falta é VPROD = 1000 — 4] 
= 959. A matriz. EXSREST é construída com os valores 
actualizados, a busca do elemento de stock mínimo não 
incidindo sobre as classes possíveis já requisitadas, neste 
caso, a classe 3 (as respectivas entradas são inicializadas a 
infinito). 


EXSD =| 450 | EXSREST =| — 1468 so 
743 — 216 — 216 

877 so — [041 

2 Do [e] 

759 — 1159 so 

591 — 2286 so 

2500 1541 —377 


Por exemplo, 


EXSREST (7,2) 


EXSD (7) — QINC (7,2) . (1000 — 41) 
2500 — 3 x 959 
= — 37 


O elemento mínimo da matriz é EXSREST [6,1 ] 
= — 2286. Na segunda iteração é requisitada a classe possi- 
vel | com o volume de produção máximo permitido pelas 
existências disponíveis dos items que a constituem. Esse 
valor é 


as0 825 759 591 2500 | 
veaspj=lLma( — — — — —, vrroD) 4 
2 tt E dg 


Como VPROD = 959, vem VCLAS [1 ] = 197. 
O volume de produção em falta é 
VPROD =959 — 197 =762. 


É escolhida a classe 2 depois de se ter actualizado o 
vector das existências disponíveis. 


EXSD =[56 546 877 2 365 0 2303 ]! 


O volume de produção permitido pela classe possível 2, é 


VCLAS [2 ]= L min ( 762 ) J = 438. 


Como não se atingiu o volume de produção total dese- 
jado e se esgotou o stock dos items do grupo, é requisitado 
o adicional numa classe possível pré-definida, por exemplo, 
a classe 3. 


O algoritmo determina os volumes de produção apresen- 
tados em VCLAS por classe possível, e requisita os items 
constituintes do grupo, nas quantidades representadas em 
QITEM. As existências disponíveis finais dos items estão 
calculadas no vector EXSD. 

VCLAS = [197 438 365]T 
QITEM = [394 1365 1971 1095 759 956 1511] ! 
EXSD = [56 000 41 0 989]! 


Note-se que 


3 
QITEM (i) = E — QINC(i;,k). VCLAS (k) 
k =1 
visto que, ; 
QINC X VCLAS = OQITEM 
(nit x nel) (nel x 1) (nit x 1) 


O algoritmo difere do caso | apenas no critério de esco- 
lha de classe possível a ser requisitada. O exemplo apresen- 
tado põe em evidência a diferença significativa que se veri- 
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fica entre os dois critérios de selecção em termos das quan- 
tidades requisitadas de material, 

É permitido alterar o critério de selecção dum grupo de 
alternativas, informando se esse grupo verifica o caso | ou 
2, por actualização da relação PRIORGRUP introduzida 
no capítulo 3 de [1]. 

De um ponto de vista prático, estes dois critérios de esco- 
lha correspondem a motivações distintas do gestor. O crité- 
rio do stock mínimo é o normalmente utilizado, pois per- 
mite eliminar os chamados “monos” em armazém — items 
em relação aos quais e por razões diversas se deixam de 
colocar encomendas. O critério de prioridades dá total fle- 
xibilidade na definição a priori de escolha. É ainda de refe- 
rir que entre iterações de processo de controlo da produ- 
ção, é possível alterar o critério de selecção. 


procedure SEL STOCKMIN; [selecção por busca do stock minimo) 


const 
nit=... (nbmero de items do grupo] 
Releases [ntmero de classes possiveis] 
infinitos... 

var 


à, j, Kel, cl, clfim: integer; 

VEALI, VrnDu, SIimih: real; 

FIM, DIF: boolean; 

EXSD, EXSREST, QITEM: array [1..mat] o real; 
VCLAS: array [1..ncl] of real; 

QINC, QINCA: array [1l..mau, I.«nci] ot real; 


function min (x, y 
[determina o mínimo eure x ey]. 


: real): real; 


begin 
liniciação de variáveis) 
for 1:=1 to nit do 
reao qtasu [i]); 
for 1:=1 to nit do 
For j:=1 tó nel do” 
begin 
read (QINC [1,9]); 
Qune L2,9):= Q1NC [4,5]; 
end; 
For 1:=1 to nit do 
ITEM [1]:=0 
for j:=1 to nel do 
tLLaS [j]J:= Infinito; 


ke): =0; 
VEALT:= VPROL; 
FIH:= False; 
ELlim:= nel: 


ijcáleuo de VCLAS] 
while not (FIM or kel=nel) do 
“begin 
kel:= kel + 1; 
SiMLH:= infinito; 
for i:=1 to nit do 
for j:=1 to nei do 
EXSREST [1,jj:= Imnilnito; 
for 1:=1 to nit do 
for j:=1 to nel do 
if not QINC [1,4]=0 
then begin 
EXSMEST [i,j]J:= EXSD [1)-GINCX [i, JJ*VFALT; 
df EXSREST [4,9] < STNLn 
tnen begin 
STMIN:=EXSREST [1,5] 
cl:=j; 
end; [then] 
end; [then] 


la viasse escolhida & cl] 
Tor i:=] to nit do 
if not QINC [i, cll=0 
then 
VCLAS veLJ:= min (VELAS [cl), [ EXSD (ij/ 
QINC [1,61]) ]; 
if not (VCLAS [cl )J=0) 


Laet Ds) lh 
1F VILAS [cl] > VFALT 
ven begin 
VELAS [clJ:= VFALT; 
FIM:= true; 


evu; (then] 
VEALT:= VEALT - VOAS [el]; 
for i:=1 to nit do 
EXSD [1]:=EXSu 11) = ViLAS [cl] * QINC [i cl]; 
end; [then 
for i:=1 to nit do 


QINCX ti,ciJj:= 0 
end; uidle] 
Tor à1:=1 to nel co 
ir VCLAS [1] = infinito 
then 
VCLas 1iJj:= O 
it not (VFALT = 0) 
then begin 
vOAS [elfim]:= VCLAS [clfim] + VFAa 
for 1:=) to nit do 
1f not (QINC [i,clfim]) = 0 
then begin 
EXSD [iJ:= EXSu La) - Venbi * QLNC [á elFfim]; 
1f EXSL (1) < 0 


then 
EXSD [il:= 0; 
end; [then] 
eia; (then] 
(requisição de material) 


for 1:;=1 to nit do 
For 4j:= to nel do 
UiTeM [1]:= GITEM [1] + VCLAS [9] * GIRC [1,9]; 


Enc, 


FIGURA 8 — Algoritmo de selecção de alternativas 
por busca de stock minimo. 


V. LOTES SIMULTÂNEOS 
23. Fases possíveis de lotes simultâneos 


Até agora, apresentou-se o estudo respeitante à monta- 
gem de um só lote. Pôs-se em relevo, porém, que vários lotes 
existem em decurso simultâneo no tempo. Nesta secção 
mostra-se que, o método de funcionamento implementado 
para um só lote, permite a existência de vários lotes em 
simultâneo, sem necessidade de complexificação. 

Em relação a um lote, começa-se por resumir O processo 
descrito, enumerando algumas regras de funcionamento. 


Facto 19: Um lote tem as seguintes regras de funcionamento 

RI: A previsão de material efectivamente fornecido é 

baseada nas existências disponíveissubsecção 11.2, 
ponto 9, facto 7). 


R2: Quando da requisição de material, as quantidades 


de material necessáriassão cativadaspara o lote 
(subsecção 11.2, ponto 10). 


Fase 1 


INFORMAT 
PRODUÇÃO 


Es ARMAZÉM 


Fase 3 


INFORMAT 
PRODUÇÃO 


EsoC) 


tg ARMAZÉM 


Fase 5 


RR = RE 
t(3) — INFORMAT 


PRODUÇÃO 
t15(2) ARMAZÉM 


Fase 2 


INFORMAT 
tg PRODUÇÃO 


ts tg ARMAZÉM 


Fase 4 


INFORMAT 
PRODUÇÃO 


t 15 (1) ARMAZÉM 


Fase 6 


t,903) INFORMAT 


PRODUÇÃO 


2/3) 


t15(5) to, (3) ARMAZÉM 


FIGURA 9 = Diferentes fases possíveis para o decorrer 
simultâneo de dois lotes. 


lote | et « lote 2 [1 
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R3: A ordem de abastecimentodas requisições de mate- 
rial não é relevante, em termos de actualização de 
dados (subsecção 111.4). 


R4: A selecção de alternativas é baseada nas existências 
disponíveis. 


Estudam-se as consequências destas regras de funciona- 
mento no decorrer de vários lotes em simultâneo. Para tal, 
veja-se a figura 9 em que se consideram os três sistemas 
INFORMAT, PRODUÇÃO «e ARMAZÉM em comunica- 
ção, no decurso da história de dois lotes simultâneos, lote 1, 


representado por MM e lote 2, representado por 1] 


Os instantes de tempo assinalados são os corresponden- 
tes aos apresentados na figura 2. 

Com o encadeamento das diferentes fases, pretende-se 
exemplificar a independência entre os diferentes lotes. 


Fase 1: O mapa de requisição colectiva respeitante ao lote | 
é abastecido pelo sistema ARMAZÉM, enquanto 
a requisição de material para o lote 2 é processada 
pelo sistema INFORMAT. 


Os fornecimentos efectivos para o lote | não são conhe- 
cidos. É a regra R2, seguida aquando da requisição do 
material para o lote |, que garante que a previsão de mate- 
rial calculada pela regra Ri fornecido para o lote 2 seja 
consistente. Por outro lado, a selecção de alternativas para 
o lote 2 é independente da situação do lote 1, pelas regras 
Rice R2. 


Fases 2, 3, 4, 5:Os dois lotes trocam as suas posições 
no tempo. 


A ordem de abastecimento no sistema ARMAZÉM é 
arbitrária, na medida em que não é descrita por nenhuma 
regra bem definida. Este facto apenas se traduz num possi- 
vel número superior de comunicação de fornecimentos 
efectivos. A análise estatística mostra porém, que este fac- 
tor não é relevante (subsecção 111.4, ponto 19). 


Fase 6:0 lote 2 termina apenas com três iterações decor- 
ridas, enquanto o lote | se encontra na 5. “iteração 
na fase de abastecimento pelo sistema ARMA- 
ZÉM., 


Embora o lote 2 tenha iniciado após o desencadear do 
lote 1, aquele termina antes deste, sem que a informação 
local de cada um dos sistemas seja por isso afectada. 

As regras RI, R2 e R4 são responsáveis pela desejada 
independência entre os diferentes lotes decorrendo em si- 
multâneo. A regra R3 liberta a informação local de cada 
sistema de relações de ordem entre diferentes lotes. As re- 
gras Ri, R2, R3e R4 traduzem a robustez do sistema do 
controlo da produção implementado. 


Facto 20: A existência de vários lotes em simultâneo ape- 
nas pode afectar o volume de dados na comuni- 
cação de fornecimentos efectivos. 


24. Comunicação entre sistemas respeitando 
vários lotes 


O sistema PRODUÇÃO pode identificar vários lotes em 
simultâneo, requerendo um mapa de requisição colectiva 
comum para os vários lotes. O processo de requisição co- 
lectiva é sequencial para cada lote, mas o mapa é comum 
aos vários lotes. 

Os fornecimentos efectivos traduzem quantidades forne- 
cidas para o somatório dos lotes. Há necessidade de frac- 
cionar estes fornecimentos pelos diferentes lotes. 


Seja QTOT, a quantidade total fornecida efectivamente, 
QREQ;, a quantidade requisitada para o lote i e QFOR;, 
a quantidade efectivamente fornecida para o lote 1. 

O sistema INFORMAT actualiza os dados de cada lote, 
calculando QFOR;, por simples divisão proporcional. 


QREQ; QFOR; 
É QREQ; QTOT 
donde, 
QREQ;. QTOT 
QFOR; = E 
í QREQ; 


Caso o item seja inserido na requisição pelo sistema 
PRODUÇÃO, a comunicação do fornecimento efectivo do 
item deve ser feita lote a lote, já que o sistema INFORMAT 
não pode calculá-lo por desconhecimento de QREQ;, para 
cada lote. 


VI. CONCLUSÃO 


Apresentou-se neste artigo a actividade do controlo da 
produção numa linha de montagem de uma unidade indus- 
trial. Descreveu-se o processamento sob a perspectiva de 
comunicação entre diferentes sistemas. Foram analisadas 
as actividades locais de cada sistema interveniente no pro- 
cesso. Os meios de comunicação e suas características fo- 
ram estudados. Discutiu-se a interacção dos processos com a 
informação e sua actualização. No capítulo 7 de [|] deta- 
lha-se o algoritmo implementado para controlo de concor- 
rência de acessos simultâneos por diferentes sistemas que 
partilham entre si os mesmos dados. 

Faz-se notar que a apresentação teve como objectivo pôr 
em evidência os pontos essenciais do processo, não-se tra- 
tando de uma descrição exaustiva. 

Para estudo do software desenvolvido, remete-se o leitor 
para o manual [3]. 
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AS VARIÁVEIS MINERALÚRGICAS 
NO PLANEAMENTO MINEIRO * 


JOSE QUINTINO ROGADO 
Professor do IST — CVRM (INIC) 


RESUMO 


Justifica-se o interesse em utilizar variáveis mine- 
ralúrgicas na avaliação, concepção e controle do pro- 
jecto mineiro. 

Apresentam-se métodos globais de optimização 
usando tais variáveis e descreve-se o método de aloca- 
ção local como forma mais conveniente de resolver 
o problema. 


ABSTRACT 


The interest of the utilization of mineralurgical 
variables for the assessment, design and control of 
the mining projects is explained. 

Several! global methods to tackle with these varia- 
bles are presented. The so called local mehod is 
develloped and it is shown to be the best way for 
solving such a problem, 


|. Tem vindo a ser reconhecida a vantagem de 
usar variáveis mineralúrgicas na avaliação de recursos 
minerais, na concepção dos respectivos projectos de 
exploração e na sua execução controlada [I, 1980, 
1984; 5, 1984; 7, 1981; 8, 1971/81]. 

Compreende-se a razão desta nova abordagem que, 
aliás, decorre da prática ainda habitual (na qual as 


principais variáveis ou atributos levados em conta são 
a mineralogia, os teores do minério tal qual, etc.), 
se se notar que, em termos monetários, a produção 
mineira é o «metal» recuperado, isto é, o metal con- 
tido in situ, diminuído das perdas inerentes à lavra 
e à concentração. 

Deixando de lado a perda resultante do método 
de lavra, considere-se apenas a perda em metal intro- 
duzida pelo processo de concentração. 

Será então (v. lista de símbolos) 

M = RmZa Q = R,Z0 

Ora, fixados que seja(m) o(s) processo(s) de bene- 
ficiação, Rm ou R, dependem de vários factores: 
Do próprio teor tal-qual, da mineralogia (composição 
mineralógica e textura) e da regulação do processo, 
quer quanto à operação de concentração, quer quanto 
à cominuição, da qual, por sua vez, depende a granu- 
lometria, uma vez fixada a textura. 


Quer dizer, o par (Rm, Z.) ou (R,, Z.) induz duas 
classes de equivalência: 

* A classe i : tipo mineralúrgico. 
* A classe J/i : ponto de funcionamento do pro- 


cesso para cada tipo mineralúr- 
gico. 


* Esta texto serviu de base à comunicação «Nota sobre Economia de Sistema Mineralúrgico» paresentada pelo autor 
ao V Encontro Nacional de Engenheiros de Minas Portugueses em Nov./84. 


Original recebido para publicação em 30/11/84. 
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Na metodologia anterior, onde as variáveis utili- 
zadas se referiam ao minério tal-qual, apareciam: 


* A classe tipo mineralógico ligada por inferên- 
cias de lógica vaga à classe 1. 


* A classe teor tal-qual da qual derivam os cri- 
térios de rejeição denominados marginalistas. 


Usando as variáveis mineralúrgicas, a abordagem 
é reconhecidamente mais directa. 

A classe ix) define antes de mais nada as partes 
exploráveis em oposição às não exploráveis. Quer 
dizer, esta classe conduz à determinação do VAB 
médio: 


VABG,))=Rm(Za,Zc/1,]) Za q— P(Zab/1,) + b(Zap/1,)) + K 
LE TA, so) 
1 E LL, m/E) 


Se o VAB,a>V, os blocos são «exploráveis», se O 


Trata-se de rejeição análoga à baseada no teor 
tal-qual, quando se utiliza o critério marginalista 
mineiro. Mas é mais rigorosa porque dá conta da 
dependência, do metal contido, de Ry, O qual, como 
já vimos, contém muita informação. 

Este critério pode ser demasiado forte, como 
aliás já acontecia com a rejeição pelo teor tal-qual. 
Tais rejeições justificam-se quando houver necessi- 
dade, por razões estruturais ou conjunturais, de tra- 
balhar com períodos de retorno muito curtos. Adop- 
tado este critério, a reserva explorável ficará definida 
e o seu complementar será eventual recurso que se 
procurará deixar in situ. 

É óbvio que, se tal não for possível, isto é, se 
houver que desmontar blocos «inframarginais», ter 
-se-á que adoptar segundo critério de rejeição, recal 
culando o VAB, deduzidas as despesas já realizadas 
por conta da reserva dita explorável [, 1975]. 

Em certos casos, como se sabe, o desmonte de 
um bloco explorável acarreta o arranque de um cone 
de influência; e o bloco só será explorável se o cone 
for «positivo» [2, 1976; 3, 1969; 6, 1979]. 

Independentemente destas considerações, há toda- 
via um ponto importante a salientar. Quando se uti- 
liza o teor tal-qual como parâmetro de selecção, a 
função critério é aditiva e a optimização consiste em 
escolher os blocos cujo teor é superior ao teor de 
rejeição; ao utilizar-se o par (My, Z.) a situação não 
é tão simples. Fixado Rm, sabe-se, como é óbvio, 
estabelecer relação de ordem entre os teores dos 
concentrados. A optimização consistirá então em 
reter os blocos com Z, superior ao teor do concen- 
trado de rejeição a R* fixado [8, 1976]. 
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Todavia, Rm ou R, podem libertar-se e o pro- 
blema põe-se na forma como se há-de construir o 
VAB optimal. 


2. Os métodos usados para resolver esta dificul- 
dade consistem em procurar pares optimais para cada 
bloco ou painel. 

Na verdade, há dois caminhos para optimizar 
tais pares. Qualquer deles exige o conhecimento das 
funções R(Z.) ou Ry(Z.) — Funções de Lavabili- 
dade. 

No método que chamaremos de alocação global, 
a função Rj(Z.)* é suposta contínua com derivada 
contínua (ambas pelo menos contínuas aos bocados) 
e conhecida num compacto limitado. 

No outro método que denominaremos de alo- 
cução local, a função Rj(Z.) é conhecida em con- 
junto finito de pontos (normalmente são mais de 
seis). Vejamos com mais algum detalhe ambos os 
métodos. 


Método de alocação global 


Como se sabe, o mercado de minérios tem sido 
dominado pelos complexos químico-metalúrgicos, do 
que resulta que, em cada caso concreto, duas situa- 
ções podem ocorrer: 


* Minérios em que o teor de rejeição do con- 
centrado se encontra abaixo do máximo da 
função Ryq; 

* Minérios em que sucede o contrário. 


No primeiro caso, há todo o interesse em vender 
concentrados com teor tão próximo quanto possível 


do teor que maximiza R,q; no segundo, devem pro- 
duzir-se concentrados com teor ligeiramente supe- 
rior ao teor de rejeição comercial. 

Assim, a alocução optimal dos minérios do pri- 
meiro tipo, faz-se maximizando o critério 

Q Ri (Z)q' (2) ie (1,...n) 

Mostra-se [9, 1972] que a alocução optimal corres- 
ponde a 


- Riq)=0 
az o 9) 


isto é, corresponde a optimizar, em cada bloco, o 
produto R,q. 

A alocução optimal para minérios do segundo 
tipo faz-se maximizando [9, 1972] 


Q, Ri 


* Ou Ro 


com a restrição 
Q Rº (Z.—-£)=0 


o que conduz ao sistema de n equações, em geral 
não linear 
dR, - 
7 PS = (1) 
R. dZi EL 
fe 


Este resultado pode escrever-se, pondo Zi — Z = 


dRi' qi 
Pp 
Atas Wi 


dr! P 


Isto é, a optimização alcança-se, alocando a cada 
bloco os pontos de funcionamento onde a elasticidade 
do rendimento em peso em relação ao teor do con- 
centrado é igual a —1. 

Claro que em ambos os casos, os algoritmos des- 
critos só serão úteis quando a função R/(Z.) tiver 
expressão analítica. 

Nas eplicações, porém, o que se conhece são 
alguns pontos dessa função, após efectuada certa 
tipologia, construída habitualmente com ensaios mi- 
neralúrgicos padrão. 

Consequentemente, haverá para cada tipo mine- 
ralúrgico uma aplicação discreta R, > Z.. Os volu- 
mes v; do tipo i a tratar com a regulação j corres- 
ponderão apenas a pontos do espaço discreto R6XZ.. 


Nestas condições, o critério será [8, 1975] 


Ri v. DE Abrs 
» e EL, eo 
com as seguintes restrições 


mi 
> Vj = D; 
J 


Ri v, (Z)—Z)=0 


Nesta forma discreta, o problema da alocução é reso- 
lúvel por programação linear [8, 1975]. 


3. Em qualquer das duas variantes (contínua e 
em geral não linear, discreta e linear) este método 
global baseia-se em tipologia prévia onde será válida 
dada relação média R(Z)? para cada um dos sub- 
jazigos tipologicamente criados. 


Aplicar este método directamente aos pequenos 
blocos de selecção conduziria à resolução de siste- 
mas de equações que poderiam ter algumas dezenas 
de milhar de equações. Com a agravante de implicar 
programação inteira. 

Em todo o caso, ficaria, deste modo, alocado a 


cada bloco o par (R,, Z.)” optimal a que corres- 
ponderia o VAB global máximo. 

Seria, portanto, possível verificar se tal VAB 
excederia certo valor mínimo fixado, Vo. 

Porém, se não se houvesse passado pela alocação 
ao nível do bloco, o resultado seria válido enquanto 
se pudesse aceitar a conjectura de que os blocos do 
mesmo tipo mineralúrgico constitutivos de um pai- 
nel mineralúrgico teriam funções de lavabilidade idên- 
ticas. E sendo estas fortemente influenciadas pela 
mineralogia, textura, granulometria e teor tal-qual, 
é evidente que estes atributos teriam obrigatoria- 
mente de figurar na matriz das contingências, evi- 
dentemente além do par (R,, Z.)“ padrão. 

Tudo isto implicaria a realização, para cada tipo, 
de krigagens e/ou cokrigagens a que se seguiriam 
amostragens compósitas para obter réplicas médias 
dos diversos tipos mneralúrgicos. 

Sobre tais amostras médias efectuar-se-iam testes 
de bancada ou piloto para construção de funções de 
lavabilidade discretas, as quais entrariam então na 
optimização global acima descrita. 

No caso de haver que rejeitar blocos, o critério 
poderia ser descartar os de menor teor tal-qual ou, 
o que provavelmente seria o mesmo, os de menor 
produto R, Z.. É bem claro que este segundo crité- 
rio seria o único correcto, se fosse possível a ele 
recorrer, o que não seria o caso ao usarem-se 
painéis mineralúrgicos formados por vários blocos. 

Porém, efectuadas as rejeições, poderia ser neces- 
sário construir nova tipologia, efectuar novas amos- 
tras médias e novos testes para se obterem novas 
funções (R,, Z.)º: as quais se utilizariam no novo 
processo de optimização (v. Fig. 1). Trata-se como 
se vê, de método iterativo, trabalhoso, com certeza 
demorado e dispendioso, que não se aconselha quando 
existe a presunção de realizar rejeição e de a fazer 
variar em face de conjuntura económica instável. 


4. Método de alocação local 


Por isso mesmo foi desenvolvido método muito 
mais apropriado que a seguir se descreve resumida- 
mente. 

Mostra-se que para constituir a função opti- 
mal Rm(S)'”, convexo superior, basta utilizar o ope- 
rador tafrach *? [8, 1974/75]. 


th Ou RL(Z). 


(2) Se no primeiro caso a fórmula de venda for linear obtém-se 


também estas equações. 

3) Ou Ri(Z). 

(44 Ou (R,,Z.)- 

(5) Ou (R,, 2.) 

(6 Ou M(C). 

71) Em homenagem ao filósofo árabe Al-Nazzam (845 dc.) 
que, para explicar o movimento, utilizou o conceito de 
salto (tafrach). 
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a 


TIPOLOGIA 
| MINERALURGICA 


KRIAGEM DE 
Ea Rm e/ou R 


COLHEITA DE 
AMOSTRAS MEDIAS 


TESTES DE BAN-| 


- Teor tq 


- Mineralogia 

- Textura 

- Resposta à cominuição/granulometria 
- Resposta a ensaio(s) padrão (RsZç) 


7 ou | Ro sZe) 


-a) max 0, Ri (7)) a! (2) 


+ E ML eseall) 


| i 1 
ou at. (Rp -q)=0 
dz 
Cc 
“, Optimiza-se, em cada tipo o pro- 


CADA OU PILOTO! 


duto Rg 
Ri (7 i 
1 =b) max Q p, 
p 
e o 
= E c/ Q ? (ei dj O 
OPTIMIZAÇÃO | K 
ou — - 7! = cte 
DO VAB dpp 
et 
dZ. 


NÃO a, > 5> SIW 


Fig. 1 - ALOCAÇÃO POR OPTIMIZAÇÃO GLOBAL 


Este operador define-se como 


AR, AR, 
E E E E 
A(RmZo)  A(R,Z.) 


E estabelece-se sobre as funções de lavabilidade 
discretas (normalmente bastam 2 a 4 pontos) de cada 
bloco, permitindo construir relação de ordem que, 
devidamente frequentada, conduz ao convexo supe- 
rior referido do conjunto de blocos e simultanea- 
mente, fornece o teor global do concentrado, que é 
assim um atributo livre. 
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O método autoriza portanto a rejeição de blocos 
quando o VAB não é satisfatório, exigindo apenas 
reordenação dos tafrach, que permitirá construção 
de novos VAB para diversos concentrados globais sem 


necessidade de utilizar qualquer combinatória. 


4. A utilização do operador tafrach, sendo muito 
mais poderosa (e portanto económica) do que os 
algoritmos de optimização global, implica, toda- 
via, a construção de alguns pontos da função Rj(Z.) 
sobre pequenas amostras e a sua subsequente esti- 
mação ao nível do bloco. 


É esta mudança de suporte que importa, portanto, 
analisar. 

Na verdade o essencial consiste em construir a 
função de libertação característica da textura de cada 
tipo mineralúrgico, o que implica análise textural de 
amostras (superfícies polidas, p. ex.) utilizando ana- 
lizador de texturas com uma janela que simule um 
grão de calibre Z que será o suporte dos histogra- 
mas de Za. 

A segunda etapa poderá ser a resolução da equa- 
ção integral 


SIF(Z, (2), dy) = F(y/2) 


Consistirá então em modelar histogramas para- 
metrizados em &, dados 


E, (Z./0) = E,(Z,): não acorre cominuição 
diferencial 
D =Y(v,v)— Y(s,s) + 0? 
e D, W Y(v, V) — Y(s, 8) 2 /s 


Esta última relação resulta de 
2 16-15) 


s/V 


E E CR ul 
o js Y6s, )— (9,0) >Y(s, s)=Y(s, 8) js 


ip Y(v, V) — Y(0,0) 


donde 


q? = vv W— Y(s,s) + o? 
o /v Y(v, V) — Y(s, 8) Gs 


Ora, exprimentalmente conhece-se, para dado tipo 
mineralúrgico, 


E, (Z,) obtido por krigagem ou krigagem dis- 
juntiva, 

—— admitindo que há variograma «in 

T(v,v; 
situ». 

TES admitindo que há variograma nas 

Ss” amostras submetidas ao analisador de 

texturas, 

q? a por estatística directa obtida no ana- 


lisador de texturas. 


A integração pode pois fazer-se por simples ajus- 
tamento de certa função de repartição, com um ou 
mais parâmetros q, que apenas dependam de 1.º e 2.º 
momentos, aos histograma sexperimentais obtidos no 
analisador de texturas, fixado q. Ou pode, obvia- 


NUMEROS DE 


mente, tentar-se resolver a equação integral. 

Seguidamente, há que construir as granulome- 
trias alargadas, dependentes, como vimos, da textura 
e do tipo e regulação da cominuição. Aqui há grande 
liberdade de métodos: Poder-se-ão usar técnicas de 
bancada ou técnicas piloto (estas de aconselhar, no- 
meadamente quando se prevê utilizar moagem autó- 
gena), etc. 

Em qualquer caso, obter-se-á, para cada tipo mi- 
neralúrgico, classes de curvas granulométricas para- 
metrizadas na composição das cargas, no tempo de 
residência, na diluição das polpas, etc. 


TIPOLOGIA 
MINERALURGIVA 


ANALISE TEXTURSL | 
HISTOGRAMAS 7,/1 | 


MODELAGEM DO HISTOGRAMA PARAMETRIZADO 
| 2 

em É, E(Z,) e Say 
(2, a(B), dy) = E(y/8) 
[8] 


> 


FUNÇÕES DE | 


REPARTIÇÃO |” LAVABILIDADE REAIS 


É e 
w 


| UTILIZAÇÃO DO TAFRACH | 


od 


CONSTRUÇÃO DO CONVEXO 
SUPERIOR M(C) 


——— — me 


REJE [ÇÃO POR 
Rota 


f 


ER 


e e 


| CALCULO DO vas OPTI ma | 


FIG. 2 - ALOCAÇÃO POR OPTIMIZAÇÃO LOCAL 
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A construção dos histogramas das granulometrias 
alargadas limita-se, como é óbvio, ao estabeleci- 
mento de combinações lineares dos histogramas mo- 
delados F(Z,, q/0), y), os quais são válidos para os 
blocos v por força de E (Z )e dq ue foram 

p ç A o) o /v q 


utilizados na determinação ou ajustamento de q. 

Vê-se, ainda, que a modelagem descrita é intei- 
ramente dispensável, quando se disponha de Fí(y/2) 
médios baseados em grandes amostras representati- 
vas e haja sido possível por krigagem alocá-lo aos 
blocos v. 

Os histogramas das granulometrias alargadas for- 
necem directamente, como é óbvio, as funções de 
lavabilidade teóricas, isto é, Fy(Z,) teóricas, por sim- 
ples introdução de critério de rejeição. 

Como passo seguinte introduzem-se as funções de 
repartição que traduzam o processo mineralúrgico e 
que, a partir das funções R(Z,) teóricas, fornecem 
as suas réplicas «reais». (Não se tratará nesta opor- 
tunidade deste importante passo do método). 

Fica-se assim em condições de aplicar o opera- 
dor tafrach que conduzirá à construção do convexo 
superior M(C) do qual se pode deduzir facilmente o 
VAB optimal fixado Z. [8, 1976]. 

E, assim sendo, as rejeições, caso haja que efec- 
tuá-las, implicarão apenas a reordenação do tafrach 
(v. Fig. 2). 

Finalmente deve salientar-se que este método es- 
trutura mudança de suporte com significado ao nível 
das variáveis mineralúrgicas. 


LISTA DOS SÍMBOLOS 


M — «metal» recuperado 

Rn — rendimento em metal ou recuperação 

Ze. — teor (ou concentração) tal-qual em metal 

O — «tonelagem» tal-gaul 

C — «tonelagem» do concentrado 

R, — rendimento em peso 

Ec — teor (ou concentração) do concentrado em 
metal 

Lib — teor (ou concentração) do concentrado em 
metal «acessoório» 

q — «valor da unidade» de metal 

VAB A.) js bruto acrescentado médio por uni- 
ade de peso tal-qual 

i — tipo mineralúrgico genérico 

n — número de tipos mineralúrgicos 

Th! — ponto de funcionamento ou processo de 
tratamento genérico, fixado o tipo i 

m/i — número de pontos de funcionamento ou 
de processos de tratamento, fixado o tipo i 

p — multa por unidade tal-qual 

b — benefício por unidade de peso tal-qual 

k — despesas de lavra, concentração, gastos ge- 
rais, despesas de capital, etc., por unidade 
de peso tal-qual 

Vo>0 — VAB médio limite económico de explora- 


bilidade 
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[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


Y; — peso do bloco ou painel do tipo minera- 
lúrgico 1 

Vij — bloco ou painel do tipo i tratado no ponto 
de funcionamento ou processo j/i 

.. — operador tafrach 

1 

A — transferência de um ponto de funciona- 


mento para outro no espaço Rx Z 
2) — calibre ou suporte 
FZ, p(9),y) — histograma cumulado (função de re- 
partição) do bloco v 

q(9) — parâmetro da função de repartição (ou da 
respectiva densidade, se existir) — escalar 
ou vectorial 

F(y/0)  — histograma médio comulado experimental, 
obtido pelo analisador de texturas fixado 
o suporte &. Ou se supõe específico da 
textura de cada tipo mineralúrgico ou por 
krigagem se faz a sua alocução a cada 
bloco v 


D. — massa do sub-depósito i 
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PROCESSOS MARKOVIANOS . 
DE RENOVAMENTO NA MODELAÇÃO 
DE SISTEMAS ESTOCASTICOS 
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RESUMO 


Este artigo consta de duas secções (parte A e 
parte B) e uma breve introdução onde se fazem 
alguns comentários sobre o método e conteúdo do 
ensino de processos estocásticos a quem tem de mo- 
delar sistemas estocásticos. 

Na parte A apresenta-se uma resenha dos resul- 
tados mais importantes da teoria dos processos Mar- 
kovianos de renovamento e uma breve discussão sobre 
as vantagens destes processos. Calcula-se a matriz 
Markoviana de renovamento, R(t), e o vector próprio 
à esquerda correspondente ao valor próprio um da 
matriz Q(+ 00), v, no caso de espaço de estados infi- 
nito numerável, especialmente, quando a matriz semi- 
-Markoviana, Oft), é tridiagonal. Neste último caso 
obtém-se uma fórmula explícita para as componentes 
de v, e mostra-se que a transformada de Laplace- 
-Stieltjes dos elementos da matriz R(t), r(s), formam 
um sistema de equações lineares recorrentes de 2.º 
ordem. No entanto, só para o «caso exponencial» 
se conseguiu encontrar a solução analítica explícita 
deste sistema. Esta solução é dada em termos de 
séries confluentes hipergeométricas. 

Na parte B tecem-se breves considerações sobre a 
literatura existente respeitante às aplicações destes 
processos à vida real; e mostram-se situações (rela- 
tivas à visibilidade de alvos que se movimentam em 
terrenos acidentados) que ilustram a maneira de usar 
a teoria dos processos Markovianos de renovamento 
na modelação de sistemas estocásticos. 


Original recebido para publicação em 7/4/85. 


24 


ABSCTRACT 


This paper begins with the presentation of some 
views on the teaching of applied stochastic processes. 

The paper is dividid into two parts. Part A pre- 
sents a summary of the Markov renewal theory as 
well as a discussion of its advantages. The main pur- 
pose of this part of the paper is to calculate the 
Markov renewal matrix, R(t), and the left eigenvector 
associated with the eigenvalue one of the matrix 
(O +co), v, in the case of infinite siate spaces, spe- 
cially, when the semi-Markov matrix, O(t), is tridia- 
gonal. In the tridiagonal case an explict formula for 
the components of the vector v is given, and the 


Laplace-Stieltjes transform of these elements of the 
matrix R(t), are expressed as a system of linear 


equations of «second order recurrence type». Explict 
analytical solution for r(s), in terms of confluent 
hypergeometric series, are obtained but only for the 
exponential case. 

In part B it is shown how to use the Markov 
renewal theory to solve a real life problem on the 
visibility of a moving target. Guidelines on the lite- 
rature of the applications of Markov renewal theory 
to real life modelling are also presented. 


INTRODUÇÃO 


Modelar sistemas estocásticos é uma arte e como 
tal aprende-se fazendo. Mas a aprendizagem é longa 
e deve ser precedida de pelo menos duas fases pre- 
liminares. A primeira consiste em aprender os con- 
ceitos básicos da teoria dos processos estocásticos 


mais usados, nomeadamente, cadeias de Markov, 
processos de renovamento, processos Markovianos, 
processos Markovianos de renovamento e Martinga- 
las. Considera-se que neste contexto o estudo da 
teoria de processos estocásticos não deve seguir uma 
abordagem baseada na teoria da medida, por esta 
ser demasiado exigente em termos matemáticos, de- 
masiado abstracta e de um modo geral dando pouca 
ênfase à interpretação intuitiva e probabilística das 
caracteríscticas fundamentais dos processos estocás- 
ticos; mas, por outro lado, também não deve ser 
feita uma aprendizagem tipo «livro de receitas» por 
não permitir uma compreensão teórica adequada que 
possibilite responder às exigências de modelação de 
problemas práticos, dada a maneira superficial e 
simplista com que aborda as questões importantes 
relativas aos processos estocásticos. O objectivo deve 
ser a compreensão dos processos estocásticos e suas 
possibilidades e limitações sob um ponto de vista in- 
tuitivo e probabilístico. 

A outra fase, que aliás vem na sequência da ante- 
rior, e a que alguns autores chamam «estudo de 
modelos», consiste em analisar e compreender as 
dificuldades que os investigadores aplicados tiveram 
de vencer, ao modelar os seus problemas concretos 
através de processos estocásticos. O estudo destas 
aplicações deve ser tomado apenas como um passo 
intermédio na aprendizagem de como fazer a mode- 
lação do sistema em estudo. É necessário examinar 
como estes utilizadores de processos estocásticos tra- 
taram «messy data» (dados sem interesse, dados de- 
sorganizados, excessivos ou insuficientes são infeliz- 
mente o que aparece com mais frequência), com- 
plexidade do sistema, sistemas para os quais a teoria 
existente é insuficiente ou mesmo inadequada, etc. 
Note-se contudo que o sistema que se está a estudar, 
de um modo geral, terá sempre particularidades es- 
peciais, por exemplo, as dependências estocásticas que 
existem no nosso problema podem diferir das depen- 
dências existentes em problemas já estudados, que 
se apresentam em tudo o mais semelhantes. Mais 
especificamente, os modelos já existentes que servi- 
ram para resolver determinados problemas concre- 
tos, não devem ser considerados como receitas ou 
descrições de como fazer relativamente ao problema 
agora a resolver. Estes modelos devem ser encarados 
como exemplos a conhecer e a estudar, através dos 
quais se podem abrir perspectivas de uma maior sen- 
sibilidade e espírito crítico, às dificuldades a enfren- 
tar e às soluções a procurar para o problema em 
estudo. 


Na sequência do que se acaba de expor, neste 
artigo, tratar-se-á a teoria dos processos estocásticos 
e suas aplicações segundo a perspectiva proposta. 

Os processos estocásticos escolhidos são os pro- 
cessos Markovianos de renovamento porque têm mos- 
trado ser úteis na modelação de sistemas complexos 


e têm a vantagem de incluir a maior parte dos pro- 
cessos estocásticos padrão: cadeias de Markov, pro- 
cessos Markovianos em geral e processos de reno- 
vamento. 

Na primeira parte (Parte A — Características e 
resultados mais importantes da teoria dos processos 
Markovianos de renovamento: vantagens e desvanta- 
gens na aplicação destes processos), começa-se por 
fazer um pequeno resumo da teoria dos processos 
Markovianos de renovamento, relacionando-os com 
as cadeias de Markov e os processos de renova- 
mento e indicando as suas propriedades mais funda- 
mentais. Apresentam-se também os prós e contras 
da teoria Markoviana de renovamento e calcula-se a 
matriz Markoviana de renovamento, R(t), e o vector 
próprio à esquerda associado ao valor próprio um da 
matriz Q(+ co), v, no caso de espaço de estados infi- 
nito numerável ,especialmente, quando a matriz 
Semi-Markoviana, Q(t), é tridiagonal. Neste último 
caso obtém-se uma fórmula explícita para as com- 
ponentes do vector v, e mostra-se que a transfor- 
mada de Laplace-Stieltjes dos elementos da matriz 
R(t), r(s), formam um sistema de equações lineares 
recorrentes de 2.º ordem. No entanto, só para o 
«caso exponencial» se consegue encontrar a solução 
analítica explícita deste sistema a qual é dada em 
termos de séries confluentes hipergeométricas. Sa- 
lienta-se ainda as limitações destes processos nas 
aplicações, devido às dificuldades de cálculo do vec- 
tor próprio v e da matriz Markoviana de renova- 
mento, R(t), no caso geral de espaço de estados infi- 
nito numerável, A Parte A foi já integralmente pu- 
blicada em Ramalhoto (1984). 


Na segunda parte (Parte B — Utilização de pro- 
cessos Markovianos de renovamento em modelação 
de sistemas estocásticos) apresentam-se situações que 
ilustram a maneira de usar a teoria dos processos 
Markovianos de renovamento (e o êxito alcançado) 
na modelação de sistemas estocásticos da vida real. 


A. CARACTERÍSTICAS E RESULTADOS MAIS 
IMPORTANTES DA TEORIA DOS PROCES 
SOS MARKOVIANOS DE RENOVAMENTO: 
VANTAGENS E DESVANTAGENS NA APLI- 
CAÇÃO DESTES PROCESSOS 


A.1 — CONSIDERAÇÕES PRELIMINARES 


Os processos Semi-Markovianos (S-M) podem ser 
vistos como processos estocásticos que generalizam 
as Cadeias de Markov (C.M.) em tempo discreto e 
contínuo permitindo que os intervalos de tempo entre 
sucessivas transições sejam variáveis aleatórias (v.a.) 
não-negativas com uma distribuição de probabilidade 
qualquer (e não necessariamente exponencial comc 
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nos processos Markovianos), a qual poderá depender 
ou não do estado de partida e do estado de chegada. 

Suponhamos, por exemplo, que temos k funções 
de distribuição F,,..., F; e que sabemos que o pro- 
cesso está no «estado jy num dado instante se a sua 
função distribuição corrente for F,. Consideremos 
ainda que a sequência de estados do processo é 
determinada por uma C.M. com uma matriz (de 
probabilidade) de transição, P=[p;], definida do 
modo seguinte: Se o processo está no estado 1 a pro- 
babilidade do estado a ser visitado imediatamente a 
seguir ser ) é dada por p;. O processo estocástico 
pontual determinado por uma sequência de interva- 
los deste tipo é chamado um processo S-M (em que 
F,=F,,; como veremos mais adiante). Muitos auto- 
res chamam a este processo pontual — Processo Mar- 
koviano de Renovamento (M.R.) — reservando o 
termo de processo S-M para o estado do processo 
como uma função do tempo. Esta é, aliás, a nomen- 
clatura seguida neste artigo. 

A importância da teoria Markoviana de Reno- 
vamento reside fundamentalmente no seu amplo do- 
mínio de aplicabilidade quer a nível teórico (no estudo 
de outros processos estocásticos mais complexos tais 
como os processos semi-regenerativos) quer a nível 
prático, como mais adiante se referirá. Situação idên- 
tica se verificava já nos processos de renova- 
mento, os quais só por si não são especialmente inte- 
ressantes mas a sua teoria (em que a equação de 
renovamento, cuja solução é a função de renova- 
mento, é capaz de explicar e fornecer todas as quan- 
tidades de interesse prático ligadas ao processo, ver 
Feller (1966)) representa um dos instrumentos mais 
valiosos dentro da teoria das probabilidades aplicadas. 


Os esforços de investigação em Processos M.R. 
orientam-se, em geral, para o desenvolvimento da teo- 
ria das equações M.R. e do modo como elas aparecem 
no estudo dos processos semi-regenerativos (i.e. pro- 
cessos que, geralmente, não são Markovianos mas 
que ainda possuem a propriedade Markoviana forte 
em certos instantes escolhidos aleatoriamente; um 
processo M.R. pode então ser «embedded» nestes 
instantes e a partir daí podem escrever-se as equa- 
ções M.R. correspondentes, cujo comportamento vem 
permitir estudar certas características de interesse 
do processo semi-regenerativo). 

As equações M.R. formam um sistema de equa- 
ções integrais que generalizam de um modo natural 
e intuitivo a equação de renovamento, dos processos 
de renovamento. A solução deste sistema baseia-se 
fundamentalmente na matriz M.R., R(t), a qual ge- 
neraliza a função de renovamento. Para detalhes ver 
Cinlar (1969, a). Este artigo «survey» é de importân- 
cia fulcral dentro da teoria dos processos M.R.; o 
método e técnicas usados neste artigo são muito 
semelhantes aos aplicados à teoria do renovamento 
em Feller (1966). 
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A teoria dos processos M.R., sob um ponto de 
vista formal, está quase completamente feita, visto 
que as questões mais relevantes estão identificadas e 
razoavelmente compreendidas. A literatura sobre es- 
tes processos estocásticos encontra-se espalhada por 
todos os periódicos de relevo ligados à chamada Ma- 
temática Aleatória (em oposição à Matemática Deter- 
minística) e suas diversas aplicações à Biologia, 
Engenharia, Medicina, Psicologia, etc. No entanto, 
só agora começa a fazer parte dos livros de texto 
sobre processos estocásticos aplicados, embora tenha 
entrado na literatura de investigação há cerca de 
30 anos e a sua aplicação a tópicos como, por exem- 
plo, Filas de Espera venha de longa data; Finch 
(1959) estudou o processo partida de um fila de 
espera MI|G|1 (processo chegada, Poisson, M, tempo 
de serviço, variável aleatória com uma distribuição 
qualquer, G, e número de servidores, 1) através de 
um processo M.R. 

Os processos Semi-Markovianos, tal como os en- 
tendemos neste artigo, foram introduzidos por Lévy 
(1954) e Smith (1955) sensivelmente na mesma altura 
mas de forma independente. Mais ou menos pela 
mesma data Takács (1954) e (1957) estudou um pro- 
cesso que sob vários pontos de vista é equivalente 
ao processo Semi-Markoviano. 

O termo processo M.R. é devido a Pyke (1961, a) 
que analisou extensivamente vários aspectos dos pro- 
cessos M.R. numa série de artigos, Pyke (1961, b), 
Pyke e Schaufele (1964) e (1966). 

É curioso notar que, até 1960) foram publicados 
apenas uma meia dúzia de artigos sobre estes pro- 
cessos mas que em 1972 tinham já sido publicados 
sobre eles mais de trezentos artigos. Esta rápida 
expansão do estudo destes processos deve-se princi- 
palmente à simplicidade da teoria, aliada ao reco- 
nhecimento do seu vasto campo de aplicações teóricas 
e práticas. 

O «Primeiro Simpósio Internacional sobre Pro- 
cessos Semi-Markovianos e suas Aplicações», reali- 
zou-se em Bruxelas, Bélgica, de 4 a 7 de Junho de 
1984; a maior parte das comunicações apresentadas 
dizem respeito a aplicações destes processos. 

Nesta secção iremos tratar de processos M.R. 
com espaço de estados, E, infinito numerável. A bi- 
bliografia apresentada até agora trata fundamental- 
mente de processos M.R. em que E é finito ou infi- 
nito numerável; no entanto, E pode ser um espaço. 
arbitrário, ver Çinlar (1969, b). Outro tipo de exten- 
são dos processos M.R. pode ser obtida fazendo com 
que estes processos representem generalizações dos 
processos com incrementos independentes estacioná- 
rios, ver Ginlar (1972) e Ezhov e Skorokhod (1969). 


A.2— DEFINIÇÃO E ENQUADRAMENTO DOS 
PROCESSOS MARKOVIANOS DE RENO- 
VAMENTO 


Seja E um conjunto finito ou infinito numerável, 
N o conjunto dos inteiros não negativos, incluindo 
o zero, e R, o conjunto dos reais não negativos. 
Sejam X, e T, as seguintes variáveis aleatórias defi- 
nidas num espaço de probabilidade (Q, F, P), 


Xº:05>E;T:05R, 


e definidas para cada neN de modo que 0= 
= 146 158% ce 


Def. 1 — Estas variáveis aleatórias X, e T, para 
n=0,1,... formam um processo estocástico Mar- 
koviano de renovamento, (X,, T,; neN), com espaço 
de estados E desde que: 


PSL TEU Te TT QU 


= PAS - | á É = I, < tX 


para todo o neN, jeE e teR 


De um modo geral, considera-se que o processo 
é homogéneo no tempo, isto é, 


PIXEL s t|X. =i]= Q,; (t) (2) 
independentemente de n. 


Def. 2-— A família de probabilidades, 


Q (0 = [Q,, (O] (3) 


para 1, jceE, teR,, dá-se o nome de matriz Semi- 
-Markoviana. (Convenciona-se que Q, (0)=0 para 
todo i, jeE). 

Note-se que para cada par (i,j)) a função: 
t— SM (t) tem todas as propriedades de uma função 
de distribuição, excepto que: 


Lj 


lim Q ()=Q, (+00) (4) 


=P [X 


n+1 


=X, = 3] ;,JeE, 


que nós designaremos por p, » não é necessariamente 
um. Por outro lado é também fácil ver que: 


p,20€ E P,,=1;LJCE, (5) 
jeE 


isto é, P=[p;] é a matriz de (probabilidade) tran- 
sição da C.M,., (X,; neN), com espaço de estados E. 


(Esta matriz e principalmente o seu vector próprio 
à esquerda correspondente ao valor próprio um, de- 
sempenha um papel de alguma importância no estudo 
de processos M.R.). 

Se E = O para algum par (1,)) então Q,, ()=0 
para todo o t e convenciona-se que (Q,, (D)/P,, = 1. 
Com esta convenção, define-se 


Q,() 
G. (1) = = 
14 (D p,. 
=P [Ti Es Tu < HXn iss I, Kn+ cs JJ; (6) 


para i, jeE e teR,. 


Ou seja para cada par (1,)) a função t— G,, (t) 
é uma função distribuição. Combinando este resul- 
tado com (1) podemos facilmente provar que as va- 
riáveis aleatórias, incrementos, T;—To, T;—T;, ... são 
condicionalmente independentes dada a cadeia de 
Markov X,, X, ..., € que a função distribuição de 
TT, depende apenas de X e X.. 

É evidente que se o espaço de estamos E tiver 
apenas um ponto, os incrementos T,— To, ... são inde- 
pendentes e identicamente distribuídos e o processo 
(T,; neN) é de renovamento. 


Um ocutro resultado importante é o seguinte: 


Seja ieE um estado fixado à priori e sejam 
Ti, Ti,... os sucessivos T, para os quais X,=i, o 
processo (Ti; neN) é de renovamento (possivelmente 
retardado). 

Portanto, a cada estado 1 corresponde um pro- 
cesso de renovamento (Ti ;neN). A sobreposição de 
todos estes processos de renovamento dá o processo 
pontual (T,;neN). O iésmo processo de reno- 


vamento contribuirá com o ponto T, se e só se 
Xn=1. Os estados destes pontos sucessivos, nomea- 
damente Xo, X1, -.., formam uma cadeia de Markov. 
O nome do processo Markoviano de renovamento 
fica portanto totalmente justificado. 

Considera-se que sup, Ti(w) = + «o para todo o 
weS. Então, para qualquer weQ e teR, existe algum 
inteiro n tal que Ti(w) <t< To. (w). Pode-se, deste 
modo, definir um processo de parâmetro (tempo) 
contínuo, (Y,; teR,), com o espaço de estados E e, 
em que: 


Y =X se É Sp (7) 


O processo (Y,.;teR,) assim definido é um pro- 
cesso semi-Markoviano com espaço de estados E e 
matriz (de transição) Semi-Markoviana Q(t) = IQ, , 
(]. | | 
Para uma melhor compreensão dos processos que 
estamos a estudar vejamos a seguinte interpretação 
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física dos processos Semi-Markovianos (Y,.;teR,): 

Pode pensar-se em Y, como sendo o estado no 
instante t de um sistema ou de uma partícula que se 
move de um estado para outro com tempos de per- 
manência aleatórios entre esses estados. O compri- 
mento de um intervalo de permanência [T,T. 
é uma variável aleatória cuja distribuição depende do 
estado X, que está a ser visitado no presente e do 
estado X | a ser visitado imediatamente a seguir. 
Os sucessivos estados visitados formam uma cadeia 
de Markov e, condicionados a essa sequência, os 
sucessivos tempos de permanência são independentes. 

A justificação do nome atribuído ao processo 
(Ye: teR,) é a seguinte: 


Se uma matriz Semi-Markoviana Q(t) é da forma: 


= Ot 
Q,0=p01-0 205; (8) 
pode provar-se que: 
PM =iYustd=PEOS = ir] (9) 


para todo o t,seR, e jeE. Mais ainda, P [Y = 
= j|Y,=1] é independente de t. 

Ou seja um processo Semi-Markoviano com Q(t) 
dado por (8) torna-se um processo (homogéneo no 
tempo) de Markov. Qualquer processo Markoviano 
regular (i.e., se para quase todos os we Q: a) t=> Y(w) 
for uma função contínua à direita e b) sup, T, (w) = 
= +oco) é automaticamente um processo Semi-Mar- 
koviano. Para detalhes ver Cinlar (1969, a). 

(Portanto, relativamente a processos Markoviano, 
a característica nova dos processos Semi-Markovia- 
nos é a da liberdade permitida na escolha das dis. 
tribuições dos tempos de permanência como já se 


tinha referido. Esta liberdade, porém, é alcançada 
à custa da propriedade Markoviana (9), a qual em 
vez de ficar assegurada para todo o t, agora fica 
apenas assegurada nos pontos de transição T). 


NOTA 1: Na definição de processos semi-Mar- 
kovianos existem algumas questões que exigem um 
tratamento cuidadoso. Uma destas questões diz res- 
peito ao facto de na definição de processos Markovia- 
nos de renovamento se permitirem transições de um 
estado para si próprio, o que faz com que tais tran- 
sições não sejam notadas por um observador de 
realizações do processo Semi-Markoviano (Y,; teR+). 
Assim quando em termos do processo Semi-Mar- 
koviano se consideram os instantes de ocorrência 
To, T;, ..., como instantes de transição e Xo Xs,... 
como os estados sucessivamente visitados é preciso 
tomar em consideração o facto de que duas ocor- 
rências sucessivas nem sempre implicam uma mudança 
de estado. 
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Uma outra questão diz respeito a problemas de 
«extinção» e de «explosão» que não referiremos 
aqui. Ver ÇCinlar (1969, a). 


A.3— MATRIZ MARKOVIANA DE RENOVA- 
MENTO 


Um processo M.R. fica, em princípio totalmente 
definido através do conhecimento da sua matriz 
Semi-Markoviana Q(t) = 1Q, (0] ad É G, (01. Em 
que a matriz IQ,, (+ 00)] = lp, ), é a matriz de 
transição da C.M. que governa a sequência de esta- 
dos sucessivos e em que G,, (t) representa a função 
distribuição do intervalo que começa no estado i e 
acaba no estado j. 


Muitas das propriedades dos processos M.R. são 
obtidas directamente das propriedades corresponden- 
tes dos processos de renovamento (lembrar que, um 
processo M.R. com um estado é um processo de 
renovamento). Por exemplo, os instantes de ocor- 
rência de um determinado estado, seja o estado ii, 
(Ti;n =0,1,2,...) formam como já se disse um 
processo de renovamento; portanto se é nestes pon- 
tos que estamos interessados, então só é necessário 
considerar a distribuição dos intervalos entre suces- 
sivos instantes de ocorrência do estado i e usar a 
teoria do renovamento. 

Em teoria do renovamento a função de renova- 
mento desempenha um papel fulcral. Em processos 
M.R. também se define uma quantidade correspon- 
dente a que se dá o nome de matriz M.R., ou matriz 
de função de renovamento e que desempenha do 


mesmo modo um papel fundamental no estudo des- 
tes processos. 


Seja N, O o número de visitas ao estado j du- 


rante o intervalo (0, t) quando é sabido que no ins- 
tante zero o processo estava no estado 1. 

Def. 3 — Seja R,, ()D=E IN, , (t)] e seja R(t) a 
matriz de ordem s (s< vc) com elementos R, (U. 
A esta matriz dá-se o nome de matriz Markoviana 
de renovamento. 

Quando s = 1, (Ny(t); t > 0), define um processo 
de renovamento e consequentemente R 1 (0) = Rt(t) 
é a função de renovamento usual deste processo de 
renovamento. 

Vejamos em seguida como é que a matriz M. R. 
está relacionada com a matriz S-M, Q(t). 


S(O=PIX =iT,<t|X, =il;i jeB; teR, 
(10) 
(para evitar trivialidades consideramos que para todo 


o 1€E, P[T;=T=.. T;=..=0|X =i]=0). 
Isto é, para n = 0: 


l,1=) 


O =I(tin= «t>b0 1 
Q (D=I(,)) ar (11) 
e para n=1,2,... 

QU) =PIX, =) Ta <tX=il 


t 
= z A dP[X=kT<u|X =). 


keE O 


rd id 


= E Í (d Q, (U)Q7, (t— u) 
ke E O 
“E WO (2 


Ou seja para que X —j e Tim <t dado que no 
instante zero o processo estava no estado i, é preciso 
que a 1.º transição para um estado k qualquer tenha 
ocorrido num dado instante u inferior a t e depois, 
começando novamente desse estado (tomando «u» 
agora como origem dos tempos) as n transições 
restantes tenham levado não mais do que «t — u» 
unidades de tempo e tenham acabado no estado j. 

Para cada estado fixo, seja j, com jcE, os ins- 
tantes T, para os quais se tem X, = j formam um 
processo de renovamento (possivelmente retardado) 
que já designamos por (Ti; = 0,1,...). 

O número de ocorrências do estado J no intervalo 


[0,t] é, obviamente dado por: 
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z o sê o 
n=0 (Ka - 1 E” < t) (13) 
(I; é a função indicatriz do acontecimento B) e por- 
tanto, 


R,,(O=EIN, (D= 


DO 
= |Xo = 1] 

n=0 Lex =iT a 4 

co 

= 2 P[X-,T;<t| X=1] 

n=0 

co car 
- 2 PYOUsijEet<o. (14) 
n=0 


Convenciona-se que o número de ocorrências num 
intervalo finito é finito, pelo que R, (t) é finito para 
qualquer 1, jeE e t< oo. Isto é, a matriz M.R., R(t), 
tem todos os elementos finitos e é dada por: 


R(D = IR, (DJ = 


5 Q (o |- = O ) 
ado ú ei hai flo) 


Def. 4 — Seja F ;€O) a função distribuição do 
«tempo de primeira passagem» do estado i ao estado 
j ou seja 


FO=PIRSHUAS] su boEls (16) 


Def. S—Seja F. (t) a função distribuição do 
«tempo de recorrência» do estado j. Ou seja, 


F (O=P[T) <t|X=il (17) 


(notar que P[Ti=0|X =3] = 1). 


Ê fácil ver através da teoria do renovamento que: 
Do 
R, (D= Ra je (18) 
e que, 


t 
R,(0= f dPM<uUX=I)Ro (tu, iz) 
0 


— far, (u) R. (t — u) 
Ê ' 


= LA ” Ru) (19) 
onde Fº (t) representa, obviamente, a n.ºsimi convo- 
lução da função distribuição F,,(.) consigo própria, 


Assim. pelo menos teoricamente, uma vez identi- 
ficada a matriz Semi-Markoviana Q(t) do processo 
deve ser possível através das fórmulas (11), (12) e (14) 
obter a matriz Markoviana de renovamento, R (t), 
e a partir desta obter as funções de distribuição dos 
«tempos de recorrência» de cada estado e dos «tem- 
pos de primeiras passagens» de um estado para outro. 

As fórmulas (14), (18) e (19) sugerem as trans- 
formadas de Laplace-Stieltjes como o método a seguir. 
As transformadas de Laplace-Stieltjes serão aqui re- 
presentadas por letra pequena, por exemplo, f ;€8) = 


=L-SIF, (0; s]= Sede, (1),s>0. 
0 

Note-se que a transformada de Laplace-Stieltjes da 
matriz Q"(t) converte-se em q"(s), isto é, a n.“ima 
potência da matriz q(s). A matriz q(s) é a matriz 
formada pelas transformadas de Laplace-Stieltjes das 
componentes da matriz Semi-Markoviana Q(t); o 
mesmo acontecendo com fi). 


Usando transformadas de Laplace-Stieltjes em (18) 
e (19) temos que, 
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00 l 
I, (s) = A E, (s) — T-L o S > Q (20) 

e 

OD=[, Or, (M;s20 1) (21) 
Ou, 

f()=1- - (22) 

ii E (s) * R 
e 

r, ,(s) | 
f,, (9) = TO 1). (23) 


Estes resultados são válidos quer E seja finito ou 
não. 

No caso da fórmula (15) a questão complica-se 
se E for infinito (numerável) pois neste caso estamos 
a trabalhar com matrizes (de funções) que são dupla- 
mente infinitas. Aplicando transformadas de Laplace- 
-Stieltjes à fórmula (15) temos que, 


Hs =I+a(s)+qX)+...;s>0 (24) 
o que implica, 

rís,=I+a(s) rs),r(s)=I+r(s)q(s), (25) 
ou seja, 

U—a(]r(s)=I,r(s) [I— a(s] =1I (26) 


Pela teoria de matrizes, se E finito então, r(s) será 
sempre dada por: 


r(s) = [1 — q(s)J" (27) 


isto é, r(s) é a matriz inversa da matriz [I — q(s)]. 
Se E é infinito numerável, então [I — a(s)]! pode não 
existir. No entanto, pode obter-se r(s) através da reso- 
lução da equação matricial, r(s) [I— q(s)] = 1, a 
qual é equivalente a um sistema infinito de equa- 
ções lineares. (Para um tratamento mais completo 
desta questão ver ÇCinlar (1969, a). 


A.4 — SINOPSE DOS RESULTADOS MAIS 
IMPORTANTES 


No caso de E infinito numerável tal como no 
caso de E finito uma vez identificada a matriz Q(t) 
e calculadas a partir dela, as componentes v, do vector 
próprio à esquerda correspondente ao valor próprio 


30 


um da matriz Q+Hco), vQ+Ho)=v, e a matriz 
Do 

M.R., Rí(t) = E Qr(t), podem obter-se todas as 

quantidades importantes associadas ao processo M.R., 

(Xw I;n = 0,1,2,...), e ao processo Semi-Marko- 

viano, (Y, t > 0), nomeadamente: 


a) tempo médio de recorrência do estado k, 7; 
Em = Á m, v), em que m, representa o 
tempo médio de permanência no estado j e é 
dado por ju — C, (D)]dt e em que, C, (1) = 


à, Q, , (b), ver Ginlar (1969, a) e 1975); 


b) distribuição limite da C.M. correspondente, 
lim P[IX =k]=vm,, ver Cinlar (1969, a) 


Tlm 9 


e 1975): 
c) média do n.º de entradas no estado k entre 
duas entradas consecutivas no estado jJ, que 
é dada por v, dividido por vw, ver Çinlar 
(1969, a); 
distribuição condicionada, P[Y,=k|X = j] 
= [1=Cit=9)] dR Oy), ver Cinlar (1969,a), 
(7.9)), se my < o pode também demonstrar-se 


o 

que lim P[Y=k | Xo=j]=(my)/( X vi mi), 
ver Cinlar (1969 a) e 1975); 

transformadas de Laplace-Stieltjes das funções 
de distribuição dos «tempos de recorrência de 
cada estado», fórmula (22): 

transformadas de Laplace-Stieltjes das funções 
de distribuição dos «tempos de primeiras pas- 


d 


e! 


e! 


C 


f 


Me? 


sagens de um estado a outro», fórmula (23). 


Existem também outras quantidades de grande 
interesse prático, nomeadamente, «tempo total de 
permanência num dado conjunto de estados durante 
um determinado intervalo de tempo finito», que 
podem ser estudadas através de um outro processo 
M.R. associado ao primeiro e cuja matriz Semi- 
-Markoviana, Q(t), está relacionada com a matriz 
Semi-Markoviana, Q (t), do processo M. R. inicial do 
seguinte modo; 


dQ(d) = é dQ(n), (28) 


ck 
onde a matriz e é dada por, 


- a E 0 0 
=P t 
O e! 0 
= t 
O O e º 


Dada a maneira simples como OK) e Q(t) estão 
relacionadas é fácil ver que uma vez calculada R(t) em 
termos de Q(t) não é difícil determinar R(t) (usando, 
geralmente, técnicas semelhantes às usadas no cálculo 
de R (t)). 


EQUAÇÕES MARKOVIANAS DE RENOVA- 
MENTO 


Sejam f,(t) e g;(t) funções não negativas, limi- 
tadas sobre intervalos finitos tais que: 


O=2(0+ = fd QuUWf(t-u),icE (29) 
kEE 


Obviamente, este sistema de equações integrais 
pode tomar a forma matricial seguinte: 


E (t) = g(t) + (Q*D (1, 30) 


onde f e g são vectores coluna cujos elementos são 
as funções f,(t) e g;(t) para 1€E. A equação matri- 
cial (30) dá-se o nome de equação Markoviana de 
renovamento (M. R.) e desempenha um papel fulcral 
na teoria dos processos M.R. e sua aplicação ao 
estudo dos processos semi-regenerativos. Muitas das 
quantidades de interesse dos processos semi-regene- 
rativos podem ser escritas sob a forma de uma equa- 
ção M.R., como veremos no exemplo apresentado 
na Parte B deste artigo. 

Quando E é constituído por um único elemento 
esta equação converte-se na «equação de renova- 
mento» dos processos de renovamento. Na maior 
parte das aplicações os processos M.R. que lhe cor- 
respondem são tais que uma vez conhecidas a matriz 
semi-Markoviana, Q(t), e o vector g(t) a equação 
M.R. (30) tem uma solução única, f (t), que é dada 
por: 


ft) = (Rºg) () (31) 


onde R (t) representa a matriz M. R. correspondente 
à matriz Q(t) conhecida. ou seja, sob a forma de 
componentes: 


É; (O) = A os du Ri(u) gu(t—u), ieE. (32) 
EE 


Embora no caso do processo de renovamento a 
solução da equação de renovamento seja única (a 
função de renovamento) no caso do processo M. R. 
não o é necessariamente. Consultar Çinlar (1969, a) 
para um estudo detalhado da equação M. R., unici- 
dade de solução, comportamento assintótico e apli- 
cação ao estudo do processo semi-regenerativo. 


A.5 — ESPAÇO DE ESTADOS INFINITO NUME- 
RÁVEL 


Quando E deixa de ser um conjunto finito e 
passa a ser infinito numerável a matriz Semi-Marko- 
viana, Q (t), passa portanto a ser uma matriz dupla- 
mente infinita. 

Neste caso mesmo o cálculo de v (VQ(+Ho) = y, 
o que exige a resolução de um sistema infinito de 
equações lineares), quando Q (+) tem todos os seus 
elementos especificados e quando a matriz é irredu- 
tível e recorrente positiva, pode não ser analitica- 
mente fácil ou mesmo determinável (dependendo tal 
facto da estrutura particular de Q(+ o), e ser neces- 
sário recorrer a métodos de aproximação por trun- 
catura. 

Para compreensão das dificuldades teóricas na 
obtenção de v e para apresentação de métodos de 
obtenção de v fundamentados em técnicas de aproxi- 
mação finita baseadas em truncaturas finitas suces- 
sivas de Q(+) e de fácil implementação em com- 
putador, ver Seneta (1981, cap. 7). 

A questão complica-se ainda mais quando se pre- 
tende obter a matriz M.R., R(t), a partir da matriz 

[8 


Semi-Markoviana, Q(t) (RC) =X, Qr(t), em que Q"(t); 
n=0,1,2,..., são definidas pelas fórmulas (11) e 
(12)). 


O método baseado nas transformadas de Laplace- 
-Stieltjes, que é sem dúvida o indicado no caso de E 
finito, no caso de E infinito numerável poderá apre- 
sentar grandes dificuldades ou não ser aplicável na 
maior parte das situações, mesmo. quando estas são 
simples como, por exemplo, quando a estrutura da 
matriz Q(t) é do tipo tridiagonal. Nesse caso as 
soluções aproximadas por truncatura que se queiram 
vir a obter devem evidentemente, ser feitas no campo 
real, isto é, utilizando a fórmula (14). 


A.5.1 — Matriz Semi-Markoviana do Tipo Tridia- 
gonal 


Nas aplicações de processos M.R. ao estudo di- 
recto ou por «embedding» de contadores estocásti- 
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cos, filas de espera, controle de qualidade, teoria de 
tráfego, processos estocásticos pontuais aplicados a 
«Neuronal Spike Trains», etc., é frequente o apare- 
cimento de um processo M.R. com um espaço de 
estados infinito numerável cuja matriz Semi-Mar- 
koviana é do tipo imigração e morte, isto é, 


MU=| 0 O O 0 
AO b) PAC) 1) 


0 PO) O p Xb 
) (O) AU) O 


(33) 


para t > O. Onde cít) e pit) são quaisquer funções, 
conhecidas, de t de i e dos parâmetros do sistema, 
que satisfazem as seguintes condições: O < et) <I 
para )j=1,2,..; O<elt)< 1 para i=0,1,2,...: 
c() +e(t)< 1 para i=1,2,...; e que tornam o 
processo M.R. regular. 

A transformada de Laplace-Stieltjes de Q(t) é 
portanto, para s > 0 dada por: 


eo ()=L-sS (P, (t);s) = 


DO 
=,S exdyt)=e(s);j=0,1,.. 
9) =L—S(o,(D;s) = 
=... =e(s);J=1,2,... - 


G,(0)=0; |k-j>1 ou j=k; j,ke (0,1,...) (34) 


Na determinação de r(s) = = q"(s) vamos usar a 
k=0 


seguinte equação chave, que se deduz facilmente da 
definição de r(s): 


r(s) [1 — q(s)] = 1 35) 


Desde que o processo M.R. seja regular esta 
equação matricial terá uma solução única, r(s). Ver 
Cinlar (1969, a). 

Notar que para i fixo a equação matricial (35) 
juntamente com (34) permite-nos estabelecer a se- 
guinte condição relativamente aos F, (9); k=0,1,...: 


3 [-c(9) = Olr, (O=1; 1=0,1,..;5>0 (36) 


k=1) 


(Convenciona-se quem, == 0). 
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A equação matricial (35), válida para qualquer 
processo M. R., aplicada à matriz q(s) definida em 
(34) conduz-nos ao seguinte sistema infinito de equa- 
ções lineares: 


ADS, (6) 5, (s) = 1 

—p, (8) Ea (s) + Fo (S)=0, (s) L2 (s) =0 
—P,.1(S) Fo,,.1(S) + To) =0,.(8) Ep (8) =0 
(37) 


Ou seja, em geral, para i fixo (1 = 0,1,...) temos o 
sistema infinito de equações lineares em r, (s): 


—P, ls) E, t+, ()— 0, ls) Es O= Je (38) 
1=0,1,...:s20 


(Convenciona-se que p.i(s)=0). 


O problema de calcular os elementos da matriz 
r(s) neste tipo de processos M. R. converte-se no pro- 
blema de resolver o sistema (38) sujeito à condição 
(36). A possibilidade de obter soluções exactas expli- 
citas, r, (9), para este sistema em termos de s, e dos 
outros parâmetros das funções c;(s) e p;/(s) depende 
inteiramente da forma particular das funções c;(s) e 
p(s) em termos de j, como é sabido da teoria das 
equações às diferenças de 2.º ordem. Por exemplo, 
se p()=i/(u+Hi+S) e c(i)=ju/Gu+HAi+s) 
(o que corresponde ao processo M.R. que governa 
toda a fila de espera do tipo MiM/wo) através da 
transformação de variáveis, t, (9) = r (9)/Qu+A+s) 


para i fixo e J=0,1,... e do uso de funções gerado- 
ras é possível resolver o sistema infinito de equações 
lineares (38) convertendo a sua solução na solução 
de uma equação diferencial de 1.º ordem de coefi- 
cientes não constantes, cuja solução é dada em ter- 
mos de séries confluentes hipergeométricas, como 
veremos em seguida. 


A FILA DE ESPERA MIM| 


Seja a matriz Semi-Markoviana, Q(t), dada por: 


Quo = Ud WHO .j=1,2,... 
PA —(u+HMt . 
Q, (D)= rea [l—e e dá J=1, 2. 
| Q,4)=0 :“li=k|>1 ou 


J=k;), k=1,2,... 


a 
Q, ()=1-e (39) 
onde À e & são números reais positivos conhecidos. 
Ou seja p()=A/(u+HA+s), jJ=0,1,... e c()=Qu)/ 
[Gu+A+s) e j=1,2,...; >0. Neste caso (38) terá a 
forma seguinte: 


A 
— K-Du+ti+s ARA Ci ni (s)— 
 G+De o 
(k+D)u+h+s E, (9) 
SO, | 449=0; s>0 (40) 


ai 


Fazendo t (s) = (1/(ku+A-+s)r, (s) e substituindo 
em (40) temos, para i fixo (1=0, , e, > 


M (D+ (ku t+, ()— 
— K+Dut, (9)=8, 5 k=0,1,... (41) 


Os coeficientes de (41) são funções lineares de k 
para k=0,1,... 

Definamos em seguida as funções geradoras y;(x) 
para 1 fixo (1=0,1,...), 


ú I, «s) 


vi) = G/x,s) = rá ds ED” = 


o 
— t (s)xk. 42 
EMO (42) 
convergentes para |x|<1. 

Note-se que, a condição (36) toma agora o se- 
guinte aspecto: 


yi(D) = G(l;s) = 1/s. (43) 


Multiplicando em (41) cada uma das «equações k» 
por x”! (k=0,1,...) somando-as conjuntamente e 
usando as funções geradoras definidas em (42) con- 
vertemos a resolução do sistema (41) na resolução da 
equação diferencial de 1.º ordem: 


“(I—)y;— [AO —x)+s] y, C)+ 
+x=0: 1=0,1,... (44) 


E 
onde y(x) = d/dx y(x) = A t AS) x*"1. A fór- 
h i = L, 

mula (43) fornece a condição inicial necessária para 
determinar completamente a solução desta equação 
diferencial. 

Estas equações diferenciais de 1.º ordem conver- 
tem-se em equações diferenciais exactas através do 
factor integrante, PR OX (q —x88/ 1, isto é, 


ue QU yo) + ([-Q+9) + 
+] y()+x) e Ma 61 (45) 


é uma equação diferencial exacta, com a condição 
inicial y;(1)=1/s. Para x>0 e s>0, 


=(A/u)(1—x) (1-x) 
NE 1 = 
PA nm CD (O) + (s/+k. 


(46) 
g((s/j+k, (s/+k+I; (1/8) (1-x); 1=0,1,... 


onde &(.,.:.) é uma série confluente hipergeomé 
trica, ver Erdelyi et al (1953). 


Pela fórmula (42) temos que: 


n! É a d” y(x) | 
nu+A+s i dx” x=0 
id NR 
LS 
= o Otra as [MP 
ç = k+l; w)] | =: 


o sb=l ds (47) 


x (e ) 
dx 
A 
onde w = a dl — x). Usando a fórmula que nos 


dá a derivada de ordem n de um produto de funções 


(Ldr[aQ). DOI] /de= = 
[dx"1), obtém-se 


o 4) [daGU] /dx! [dib(x)] / 


E na 
o nuti+s Jd NO dv E) 
O E Dam 


di dC +R+I SW) | : hiMd 
e W = a . (=) “ 
(48) 


Notando que, 


d(we) k(Kk—1... [Kk—(G—Dlw'!,ifk=), JH... 


dw! 0 (49) 
, caso contrário; 


e usando ainda o resultado em Erdélyi [53; p. 254, 
(10)], para derivadas de séries confluentes hipergeo- 
métricas, temos finalmente, 


! min(n,k 


) 
— Di(t 
a ) DO 


“[a) Jetta dt 


| sq É gia Alt: 
é (n Hai; Eh 


i 
r ()=(nu+i+s) 3 (— DC) 
k=0 


ne (0,1); 3 =>0. (50) 


A.6 — CONCLUSÕES 


1. Sob o ponto de vista probabilístico e das apli- 
cações correntes é razoável convencionar, que na 
matriz M.R., Q(t), definida pela fórmula (33) se 
tenha (HO) =m=1 e 
p(+o+e(+o)=p+o,=1, para i=1,2,.... 


Esta convenção adicional, juntamente com a estru- 
tura particularmente simples da matriz Q(t) faz com 
que as coordenadas v;; k = 0,1,... do vector próprio 
à esquerda correspondente ao valor próprio um da 


matriz (+ =), (solução do sistema X v Q, (+= 
i=0 : 


v; ]=0,1,...) sejam sempre determinadas e dadas 
por: 


vw =1 


Pi. oPk-s 
= =... + = aa 1 
"= ego SS bê di 


2. No caso do cálculo da matriz M.R., R(t), a 
situação já não é tão agradável, pois embora se saiba 
que ela existe, só nos foi possível resolver analitica- 
mente a equação matricial (35) no caso em que Q(t) 
é definida por (39). A fórmula (50) dá a solução 
encontrada que aliás está longe de ser simples. 


3. Podemos portanto concluir que para processos 
M.R., de espaço de estados infinito numerável, é 
sempre possível calcular exacta ou aproximadamente 
as quantidades associadas ao processo, relativas ao 
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comportamento do estado estacionário do processo, 
isto é, que sejam do tipo do tempo médio de recor- 
rência de qualquer estado, 74; da distribuição limite 
da C.M. correspondente; da média do n.º de entradas 
num dado estado entre duas entradas consecutivas 
num outro estado, e do Ea PIY,=k | X=il. 


Visto que, estas quantidades são definidas essencial- 
mente à custa das componentes v, do vector próprio 
à esquerda correspondente ao valor próprio um de 
Q(+ 0), e são independentes dos elementos R, (b da 
matriz R(t). O que aliás, não é o caso relativamente 
às distribuições condicionadas, P[Y,=k | X, = ]l; 
transformadas de Laplace-Stieltjes das funções da 
distribuição dos «tempos de recorrência de cada 
estado» e dos «tempos das primeiras passagens de 
um estado a outro», que como vimos, fórmulas (22) 
e (23), dependem dos elementos r, (8). 


4. Em conclusão há ainda muito a fazer na área 
da aplicação da teoria às situações práticas, uma vez 
que o cálculo da maior parte das características 
importantes destes processos, (exceptuando as carac- 
terísticas relativas ao comportamento estacionário) se 
baseiam na matriz M.R., R(t), a qual é muito difícil 
de calcular analiticamente, a não ser no caso de uma 
estrutura particularmente simples da matriz semi- 
-Markoviana, Q(t). Os caminhos propostos para ultra- 
passar estas dificuldades consistem na procura de 
algoritmos computacionais e no desenvolvimento de 
teorias de aproximação para o cálculo de algumas 
destas características. Evidentemente, que estas duas 
alternativas estão fortemente dependentes das aplica- 
ções em vista. É aqui, na ligação da teoria à prática 
que se encontram as áreas priveligiadas de investi- 
gação actual em processos M.R., visto que como já 
se referiu atrás a teoria dos processos M.R., sob um 
ponto de vista formal está praticamente concluída. 


PARTE B— UTILIZAÇÃO DA TEORIA DE 
PROCESSOS MARKOVIANOS DE 
RENOVAMENTO NA MODELA- 
ÇÃO DE SISTEMAS ESTOCÁS- 
TICOS 


Os processos M.R. têm sido aplicados com sucesso 
na modelação de muitas situações práticas, nomeada- 
mente: 


— Filas de espera, Cheong, de Smit e Teugels 
(1974), Disney (1975) apresentam uma vasta biblio- 
grafia sobre este assunto se bem que não totalmente 
actualizada, visto que, nestes últimos dez anos muito 
se tem publicado sobre esta matéria; por exemplo, 
Ramalhoto e Girmes (1977), Ramalhoto (1981) estu- 
dam, numa perspectiva unificadora, filas de espera de 
capacidade infinita de serviço, através de processos 
M.R. de uma forma exacta e, também, por aproxi- 


mações consoante a natureza probalística do «input» 
e do serviço da fila; Cole (1979) apresenta um mé- 
todo interessante de amostragem, numa aplicação de 
um processo M.R. multidimensional ao estudo de filas 
de espera em aeroportos; Powers e Lackey (1975), 
Purdue (1975), estudam tipos específicos de filas de 
espera, surgidas no estudo de sistemas biológicos, em 
termos de processos M.R. 

— Transportes, por exemplo, Daganzo (1975) 
estuda problemas de tráfego rodoviário através de 
processos M.R. com três estados e concluí que o 
seu modelo parece produzir resultados mais realistas 
do que os modelos existentes até àquela altura. 

— Redes de ambulâncias, por exemplo, Hall (1969) 
apresenta um longo relatório sobre o estudo dos 
serviços de emergência da polícia de Detroit, U.S.A., 
onde trata as várias fases do problema, isto é, recolha 
de dados, elaboração de um primeiro modelo proba- 
bilístico grosseiro, estudo estatístico das caracterís- 
ticas mais importantes do sistema (de modo a fazer 
por exemplo, um ajustamento razoável das distribui- 
ções de probabilidade que é necessário postular no 
modelo probabilístico escolhido, etc.); o sistema global 
do serviço de emergência foi estudado como sendo 
uma fila de espera com um processo chegada Marko- 
viano de renovamento. Este modelo de fila de espera 
é consistente com os dados observados, fornece uma 
estrutura mais simples que muitos outros modelos 
probabilísticos previamente explorados para calcular 
as características de interesse prático do sistema, e 
fornece ainda uma melhor compreensão dos meca- 
nismos do sistema na sua totalidade. Este relatório, 
se bem que não muito recente é no entanto um 
valioso trabalho de inter-acção entre teoria de filas 
de espera, teoria de processos M.R., modelação esto- 
cástica de problemas da vida real e estatística. É inte- 
ressante notar que foi o estudo teórico e puramente 
analítico e abstracto sobre este tipo de filas de espera 
com processos chegada M. R. desenvolvido pela Prof. 
R. L. Disney, cinco anos antes do início deste pro- 
jecto, que permitiu modelar a situação prática con- 
siderada. 

— Situações referentes a gestão de pessoal e socio- 
logia por exemplo, Bartholomeu (1982). 

— Estudos de fiabilidade de equipamentos, por 
exemplo, Di Marco (1972). 

— Situações referentes às ciências computacionais, 
à psicologia e às ciências biológicas, por exemplo, 
Sim (1976). 

— Em «design» de estruturas de cimento armado 
para reactores nucleares, ÇCinlar, Bazant e Osman 
(1977). 

A literatura sobre as aplicações à vida real dos 
processos M. R. é vasta e ainda mais dispersa e menos 
acessível que a literatura sobre a teoria destes 
processos. Na verdade muitas destas aplicações confi- 


nam-se a relatórios privados e algumas vezes confi- 
denciais de empresas industriais e comerciais, hospi- 
tais, instituições de investigação, departamentos de 
defesa e do Governo em geral, assim como, em teses 
de mestrado e doutoramento em departamentos de 
Matemática Aplicada, Estatística, Investigação Ope- 
racional, Engenharia Industrial, Civil, Mecânica, 
Electrotécnica, Química, Psicologia, Sociologia, Ciên- 
cias Médicas, Biologia, Gestão, «Management Scien- 
ce», etc. 

Apresenta-se em seguida um exemplo, de uma 
situação da vida real, que se espera ilustre razoavel- 
mente a maneira de aplicar a teoria dos processos 
M.R. a situação concretas. Escolheu-se este exemplo 
por ser suficientemente pequeno para poder ser resol- 
vido aqui e permitir ilustrar o uso de muitos dos 
conceitos e dos resultados sobre teoria dos processos 
M.R. apresentados na Parte A, deste artigo. Sob o 
ponto de vista prático este exemplo se bem que trivial, 
foi no entanto usado como parte de um modelo de 
combate desenvolvido em Miller, White e Thompson 
(1979) para fornecer ensinamento sobre o uso de 
armas contra alvos móveis em terrenos agrestes. 


Exemplo sobre Visibilidade (Miller, White e 
Thompson (1979)). 


Considera-se que o alvo se movimenta em ter- 
reno acidentado. A ideia básica é a de que o alvo 
pode estar apenas em um de dois estados: invisível 
ou visível, isto é, espaço de resultados E=(0,1). 
O objectivo é estudar a visibilidade do alvo. Será o 
alvo visível no instante t? No caso afirmativo, por 
quanto tempo? Aliás uma das questões práticas ful- 
crais consiste em estudar as armas existentes ou a 
desenvolver, com o objectivo de comparar o tempo 
necessário de pontaria e chegada ao alvo do pro- 
jéctil, com o tempo de visibilidade do alvo. 

Seja, 


cai * LR od à ki me 
O, se o alvo não é visível na n transição 


Xa=— 


a 
- transição 


|, se o alvo é visível na n 


La 


Dadas as características do sistema não é possí- 
vel distinguir se um alvo visível muda de estado mas 
continua visível ou sendo invisível muda de estado e 
continua invisível e por isso assume-se que: 


O, se i=) 
=X =i] = (52) 
1, se ii, 
ou seja, se num dado instante o alvo está no estado 1, 
no instante seguinte em que muda de estado tem de 


ir para o estado j5Ai. Define-se T, como sendo o 
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“* ésima 
instante em que o alvo muda de estado pela n —— 
vez e considera-se que: 


PIT. T<UX. =X =]=6,(0 (53) 


Se o terreno tem esconderijos é razoável esperar 
que em média o alvo permaneça escondido durante 
longos períodos de tempo. Se pelo contrário o ter- 
reno é escampado então esperar-se-á que em média 
o alvo fique visível por longos períodos de tempo. 
Uma ou outra situação reflecrir-se-á no valor de 
Golo. 

Toma-se o processo base de modelação deste sis- 
tema estocástico como sendo o processo M.R. 
(X, Ta; neN) com E = (0,1) e Q(t) = [QD] = 
= [PIX,=LE TIS) =] = IPI, 
A = PII. 


+1 e 
= 1, SUA db Sado! = 
ésimo 


+1 


[1 — 1) G,()] onde 1 é o (1,)) — elemento da 
matriz identidade e convenciona-se que T;=0 e 
P[Xo= 0] =1. Isto é, 


Q(t) = 0 Gol) 7) 
| (54) 
Gt) 0 | 


(um processo M.R. com uma matriz semi-Marko- 
viana tão simplificada quanto esta goza de proprie- 
dades especiais e chama-se um processo de renova- 
mento alternado; note-se que sabendo em que estado 
se está, o estado seguinte fica completamente deter- 
minado). 

Por hipótese, o processo é considerado homogé- 
neo, isto é: 


PO =P =) T<tX =i 


PR iTo- LEtixme=en. 65 
Ou seja, pelas fórmulas (54) e (55) 


t 
Qu (O = JE Í di Qu) Qult—u) 


| 


f à Qulu) Qut—u) = 
O 


II 


(Goa * Go)(O, 
e claramente, 
QUO (1) = 0 
Portanto se o processo começa no estado “0”, as 
subsequentes entradas no estado” 0” (invisibilidade) 


formam uma sequência de v.a. (S,; neN) e (Saw Sa; 
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neN) é uma sequência de v.a. independentes e iden- 
ticamente distribuídas com 


PIS 


SS SU=(6,"* OW: (56) 
Ou seja, as sucessivas entradas no estado “0” se o 
processo começa no estado “0” é um processo de 
renovamento com função distribuição dada por (56). 

Para prosseguir o estudo da visibilidade do alvo 
é necessário introduzir outras famílias de variáveis 
aleatórias (processos estocásticos semi-regenerativos) 
ligadas ao processo M.R. inicial. 


Sejam, 
Y=A sell <t<T.rcomteR - (57) 


em que Y, diz-nos simplesmente se o alvo está ou 
não visível no instante t, e 


L=Tart' T<t<To. comteR, (58) 


em que Z, diz-nos por quanto tempo o alvo perma- 
necerá no estado que ocupa no instante t. 


Seja, 


P[Y=j Z>y|X=1]=H0,y,0. 69) 


Notar que H,(J,y,t) contém toda a informação de 
visibilidade do alvo. Por exemplo, para j=1 e 1=0 
ela representa a probabilidade que o alvo seja visí- 
vel no instante t e que permaneça visível por mais 
do que y unidades de tempo dado que inicialmente 
o alvo estava invisível. Vamos ver como obter esta 
probabilidade em termos dos elementos da matriz 
semi-Markoviana, Q(t), definida por (54) e da sua 
correspondente matriz Markoviana de renovamento, 
R(t). Para tal, e porque estamos interessados na proba- 
bilidade de ocorrência de acontecimentos num instante 
futuro t>0, vamos recorrer a uma partição baseada 
no instante de ocorrência da 1.º transição depois do 


instante zero, T;, de modo a essa probabilidade poder 
ser obtida em termos de dois acontecimentos mutua- 
mente exclusivos: t<T, e t>T, Com efeito, dado 
que X,=0, quando T,>t teremos, 


Y, =), com probabilidade 1, se j=0 e 
com probabilidade 0, se j0. 


Além disto, T>t e Z>y, se e só se, não tiver 
havido nenhuma mudança de estado antes de € (e 
portanto no instante t, Y,=0), e não vier a haver 
nenhuma mudança de estado durante as seguintes y 
unidades de tempo depois de t. Em geral, 


= cs TT Ti —— 


P [Z.>y, Ti>t | X=i]l = P[T>tty | Xo=i) = 
eg Sl ” 


e portanto, 
P [Y,=j, Z>Dy, Ti>t | Xo=il = 


=Iy [> Qu(t+y). 


Recorrendo agora à partição em T, temos: 
H, (3, LP t) = P [Y.=], Z>y, T,>t | Xo=1] 
+ P[Y,=j, Z>y, Ti<t| Xo=i 


= 1 [1 — sa Qu (t+ y)] + 


t 
+ X Sd P(X-k T<u|X=1) 
keE O 


PY, =) Z,>7 | X-H 


t 
E z s dy Qu(u) Hi), yt—u), 
keE O 


para i€cE. (60) 


A equação (60) é uma «equação Markoviana de 
renovamento», e portanto, 


t 
H, 0, J, t) = E] dy R;(u) ji-= 
- i(tty— 61 
vd On (E y—u)] (61) 


onde Rit) é o (i,j) ésimo termo da matriz Marko- 
viana de renovamento do processo M.R. subjacente 
e definido através da matriz semi-Markoviana dada 
por (54). Como vimos na Parte A deste artigo em 
geral esta matriz R(t) não é fácil de calcular. Porém 
neste nosso exemplo, dada a forma extremamente 
simples da matriz Q(t), não vai ser difícil calcular 
R(t). 

Com efeito, a transformada de Laplace-Stieltjes 
de Q(t), a(s), é dada por, q(s) = O gu(s, | em 

edi) 0 ] 

que g;; é a transformada de Laplace-Stieltjes de G;j(t). 
Como sabemos E=(0,1) e portanto finito e pela 
teoria dos processos M.R,, 


r(s) = [I—q(s)]" = 1 Bo(s) 
1I— gols) gls) 1—ga(s) Eros) | 
Eis) | 
I—go(s) gls) 1— gols) Bis) 
(62) 


(Notar que se as funções distribuição G;j(t) não tive- 
rem transformadas de Laplace-Stieltjes, não se 
poderá obter r(s)). 

Se r(s) existir, pelo menos formalmente, podemos 
inverter esta transformada e obter R(t). E como 
sabemos, pela Parte A deste artigo uma vez cal- 
culada R(t) está-se apto a resolver quaisquer tipos 
de questões relacionadas com o processo M.R. em 
causa. 

Consideraremos agora que Ga()=Ga)=l-—es!, 
(esta hipótese é pouco realista, mas é aqui conside- 
rada para simplificar os cálculos. Num problema prá- 
tico era preciso fazer uma análise estatística cui- 
dada, de modo a permitir escolher para distribuições 
Go(t) e Gy(t) aquelas distribuições que «melhor se 
ajustassem» aos dados observados, através do uso da 
teoria de testes de hipóteses e estimação), então 


H 
Bo(s)=giks)= RE e T(s) = 
E 
l | : u+s | 

|—pZ/(u+s)? H 1 
L+s 
o que não é difícil de inverter. 
| » t=0 
Ro(t) =Ru(t) = ut 1 

a LO) t>0 


t 
Rol) = Rad) = E + (mem) DO. 


l 


(Notar que R;í(t) terá sempre um salto na origem de 
tamanho 1, visto que se definiu esta função no inter- 
valo fechado [0,t] e se considerou a ocorrência de 
um salto no instante t=0 (T,=0). Portanto o número 
esperado de visitas ao estado i é sempre pelo menos 
um se X,=i). 

Através das fórmulas (61) e (63) e depois de algu- 
mas manipulações algébricas simples, temos: 


Ho(0,y,t) = S e + OEA, (64) 


HXLy) = em — Deo, (69) 
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Ou seja, a probabilidade de que o alvo que está 
inicialmente invisível, fique invisível no instante t e 
permaneça invisível durante pelo menos mais y uni- 
dades de tempo é dada por (64). A fórmula (65) dá 
portanto a probabilidade de que o alvo que inicial- 
mente está invisível fique visível no instante t e 
permaneça visível durante pelo menos mais y unida- 
des de tempo. 


H(0,y,) + H(L,y,D=PIZ>y|Xo=0]=e” (66) 


Esta fórmula representa a probabilidade de que o 
intervalo de tempo a partir do instante t até à mu- 
dança de estado seguinte seja maior do que y unidades 
de tempo, qualquer que seja esse estado seguinte que 
venha a ocupar. Visto que os tempos de ocpuação 
dos estados são variáveis aleatórias independentes e 
idênticamente distribuídas, (66) é uma consequência 
da «falta de memória» das exponenciais. 


Se em (59) fizermos y—0, teremos, 


ou seja, as probabilidades dos estados (dependentes do 
tempo) para o processo semi-Markoviano (Y,, teR,). 

Se em (64) e (65) fizermos y>0 e t=>, obtemos 
as probabilidades estacionárias dos estados para 


(Y, teR.), as quais são obviamente E porque nós 
consideramos Gal) = Guít). 

A probabilidade do alvo permanecer visível du- 
rante mais y unidades de tempo dado que neste 
momento está visível, é dado por: P[Z,>y|Y,=1] = 
=euy, t>0. Visto que: 


HXLy0)=P[V=1, Z>y|K=0]= 
1 


l 
=>— gu—— gruteuy 
2 


2 
H(L,O,)=P[Y,=1/X,=0] = — ut), 


Portanto, P[lZ>y | Y=1, X=0]= 


- PIZ>y, Y,=1]X,=0] di e 
=> Plv=iDo=o Ye como P[X,=0]=1, 


segue-se o resultado pretendido. A simplicidade 
de todas estas fórmulas resulta de se ter escolhido 
Goa(D=Giu(t)=1—e“.Neste caso (Y, teR,) é um 
processo Markoviano e (X,, Ta; neN) é um processo 
de renovamento alternado e não era necessário usar 
as ferramentas mais sofisticadas dos processos Mar- 
kovianos de renovamento, o que aliás fazemos aqui 
apenas com o intuito de ilustrar a teoria apresentada 
na Parte A deste artigo. 
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RECURSOS HUMANOS NACIONAIS 


E DESENVOLVIMENTOS 
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Director do Instituto de Energia do LNETI 


RESUMO 


O estudo apresentado consta de três partes. Na 
primeira, refere-se, sinteticamente, a situação educa- 
tiva portuguesa e comentam-se as projecções sobre 
necessidade de recursos humanos no pais até 1992, 
feitas em trabalho realizado pelo Ministério da Edu- 
cação com a colaboração do Banco Mundial. Na 
segunda parte, apresenta-se o panorama nacional nos 
domínios da LD&D comentando-se os valores dispo- 
níveis. Fecha o estudo um capítulo de considerações 
finais. Aí, salientando-se que o país não dispõe de 
recursos relevantes de matérias-primas e energéticas, 
conclui-se que se deverá investir fortemente nos seus 
recursos humanos, beneficiando da posição geográ- 
fica e das relações sócio-económicas tecidas ao longo 
da História. Consequentemente, o fortalecimento dos 
recursos humanos nacionais é a prioridade das prio- 
ridades. 


| — A SITUAÇÃO NACIONAL E O DESENVOL- 
VIMENTO 


Sendo Portugal um país sem grandes recursos de 
matérias-primas e energéticas, dever-se-á investir 
fortemente nos seus recursos humanos, beneficiando 
da posição geográfica do país e de relações sócio- 
-económicas tecidas ao longo da História. 


O primeiro investimento nos recursos humanos no 
nosso país terá, necessariamente, de ser feito na 
Educação. Todavia este esforço deverá ser norteado 
pela ideia base de que o sistema educativo tem de 
ser aberto, global, coerente, contínuo e flexível não 
permitindce que o Homem adquira posicionamento 
fixo na sociedade, resultante da educação recebida 
apenas em dada época da sua vida. Trata-se de insti- 
tuir a educação permanente que consiste em um 
sistema contínuo de formação que procura responder 
às exigências da realização pessoal e social de cada 
pessoa e da sociedade na sua evolução temporo-espa- 
cial educacional, dos jovens aos adultos, feita equili- 
brada e conjuntamente ao longo de toda a vida e por 
toda a nação. 


É preciso não só formar e informar para o 
primeiro emprego, evitando a inadequação entre o 
que se aprende e o que a vida real exige, como 
ainda formar de modo dinâmico de maneira a que, 
a todo o tempo, a pessoa possa reagir, adaptando-se 
a novas situações. Por um lado, há que obter boa 
correlação entre a escola e o mundo do trabalho, 
por outro, há que ter em conta a evolução tecnológica 
que pressupõe reciclagens e o ajuste aos novos 
saberes. 


Tenha-se em atenção que a educação é o processo 
que permite a integração da pessoa na colectividade 


* Comunicação apresentada ao Congresso da Ordem dos Engenheiros de 1985, em Coimbra. 


original recebido para publicação em 11/4/85 


40 


à medida que se desenvolve, favorecendo o desabro- 
char das suas faculdades. 


As modificações profundas da sociedade actual 
não são apenas quantitativas, mas também qualita- 
tivas. A cibernética tende a substituir-se a certas for- 
mas de actividade intelectual, a telemática permite 
ao Homem colocar-se à escala da velocidade plane- 
tária. A mutação da sociedade vai no sentido do 
terciário. O primário torna-se residual, o mínimo para 
garantir a subsistência da nação. A população migra 
para o secundário e, cada vez mais, deste para os 
serviços, concentrando-se na indústria da informação 
— meios de comunicação social, informática, tele- 
mática, ensino. 

A cultura deixou de ser veiculada principalmente 
pelo livro para o ser pela imagem, pelos educadores 
sem rosto como a televisão, o cinema, a publicidade, 
os discos, o telefone, a banda desenhada, etc. À pro- 
liferação dos meios de transmissão da cultura veio 
juntar-se a enorme massa de informação disponível 
o que exige o contínuo reciclar dos nossos conheci- 
mentos. 

É imperiosa a inter-relação entre a Educação, O 
Trabalho e a Produção. Tem de se acelerar o desen- 
volvimento e criar emprego. O desenvolvimento cor- 
responde a um estádio organizativo da sociedade 
caracterizado pela existência de acesso real de toda 
a sua população às oportunidades do nosso tempo, 
destacando o emprego e os bens e serviços associados 
às necessidades básicas. 

Neste contexto interessa analisar, ainda que bre- 
vemente, a panorâmica da situação educativa portu- 
guesa. Esta sintetiza-se no seguinte quadro de situa- 
ções: 


a) elevada taxa de analfabetismo (=23%); 

b) baixa taxa de escolarização do ensino pré- 
-escolar (= 12%); 

c) insuficiente cobertura da escolaridade obriga- 
tória (=70% da população escolarizável); 

d) baixa duração média da escolaridade obriga- 
tória (6 anos); 


e) baixa produtividade geral do ensino com taxas 
de repetência elevadas; 

f) disparidades regionais acentuadas, reveladas 
quer pelas taxas da repetência, quer pelas de 
abandono; 

g) rede escolar com muitas insuficiências; 

h) deficiente equipamento de muitos estabeleci- 
mentos de ensino; 

1) número escasso de professores profissionaliza- 
dos; 

]) recursos financeiros insuficientes postos à dis- 
posição do sector; 

|) inexistência generalizada de cursos profissio- 
nais médios; 


m) grandes assimetrias quanto a frequência e 
proveniência sócio-económica da juventude 
universitária; 

n) discreto financiamento da investigação cien- 
tífica e tecnológica (=0,35% do PIB); 

o) escasso potencial humano científico-tecnoló- 
gico do país (cerca de 2,7 investigadores por 
10 000 habitantes); 

p) fraca participação do sector empresarial nas 
actividades de I,D&D. 


Convirá ter presente que cerca de 45% dos jovens 
não completa o 9.º ano de escolaridade e que cerca 
de 67% dos jovens não prossegue os seus estudos 
para além do 9.º ano. Apenas 8 a 10% dos jovens 
segue estudos superiores. 

Temos uma situação pouco invejável em todos os 
níveis de ensino, desde o primário, pelo secundário 
até ao terciário. Pode discutir-se por onde começar, 
se pelo primário, base de todo o edifício, se pelo 
terciário procurando impulsionar rapidamente o 
desenvolvimento técnico-científico do país, suporte de 
políticas industrial, agrícola, sanitária, etc. coerentes, 
se pelo secundário terminal que apresenta lacunas 
graves entre nós. 

Creio que as três frentes têm de ser cuidadas 
simultaneamente já que a delicadeza da situação e 
a eminência de roturas graves a isso obrigam. 

O Ministério da Educação publicou em 1981 um 
estudo sobre a «estratégia de desenvolvimento dos 
recursos humanos». Este estudo foi elaborado em 
colaboração com o Banco Mundial. Procurou-se 
determinar as necessidades da mão-de-obra e a evolu- 
ção do mercado de trabalho no nosso país para o 
período 1980-1992. 

Usou-se o chamado «Modelo integrado de previsão 
das necessidades de mão-de-obra». Trata-se de um 
modelo computorizado que usa indicadores dos siste- 
mas económicos e educacional bem como da oferta 
e da procura de mão-de-obra. 


O bloco económico foi projectado com base no 
PIB por sectores de actividades e na produtividade 
média da mão-de-obra no ano de base. O bloco 
educacional faz o acompanhamento de cerca de 1200 
grupos de alunos (100 cursos X 2 sexos X 6 grupos 
etários) ao longo do sistema escolar. A simulação 
do mercado de trabalho baseia-se nos «outputs» de 
blocos anteriores que dão os fluxos que afectam 
as disponibilidades de mão-de-obra para o ano base. 

O modelo usado baseia-se na matriz do emprego 
por sectores de actividade e por grupos profissionais 
e na afectação dos alunos provenientes do sistema 
educativo a essa matriz. 

As projecções realizadas tiveram a tripla preo- 
cupação seguinte: 


— manter um ritmo de crescimento do PIB a 
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QUADRO 1 


Comparação entre emprego potencial e emprego efectivo em 1992 


Excesso (+) 
Défice (—) 
(1) — (2) 


Emprego 
NÍVEIS DE PROFISSÕES efectivo (1) 


Emprego 
potencial (2) 


Pessoal superior científico e técnico 


Engenheiros e Arquitectos 34,6 
Direito, Contabilidade e Economistas 36,0 
Professores 109,9 
Quitro pessoal superior 3241 
Quadros médios 85,3 
Empregados altamente qualificados 413,1 
Trabalhadores manuais altamente qualificados 228,4 
Empregados qualificados. Nível intermédio 190,7 
Operários qualificados. Nível intermédio 252,1 
Empregados semi-qualificados e trabalhadores agrícolas 161,1 
Operários semi-qualificados 788,5 
Não qualificados 424,3 


* Possíveis excedentes de mão-de-obra. 
Fonte: REIS A. (1981), 


QUADRO II 


Ensino Politécnico (valores médios anuais) 


coutputs» para satisfazer 
as necessidades de 


coutputs em 1992 
com manutenção 


toutputs» em 


EasA da política actual mão-de-obra em 1992 
Contabilidade e Administração 348 650 a 750 
Administração Pública ss 1000 a 2 000 
Agronomia E 600 a 700 
Eng* Alimentar 74 350 a 450 
Eng." Têxtil 8 300 a 400 
Artes Gráficas tm 100 a 150 
Eng." Mecânica 218 800 a 1000 
Eng." Electrónica e Telecomunicações 270 700 a 900 
Eng. Civil 188 600 a 900 
| Engº Química 300 a 400 
Eng* de Pesca —. 50 a 100 
Artes Decorativas e Design — 200 a 350 
Línguas e Turismo 301 600 a 700 


Fonte: MEU (1981) Bol. Inf. n.º 14. 
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níveis comportáveis com a necessidade de dimi- 
nuir a lacuna que nos separa das capitações 
europeias. 

— modernizar a estrutura profissional portuguesa 
tendo em vista colmatar, na presente década, 
cerca de metade da distância que nos separa 
da média europeia. 

— admitir que a escolaridade obrigatória passará, 
até ao fim dos anos oitenta, de 6 para os 9 anos. 


Embora as hipóteses de partida se tenham reve- 
lado optimistas, os valores encontrados e que se apre- 
sentam são os únicos disponíveis e funcionam como 
indicadores desejáveis. 

No (Quadro I é feita a comparação entre o 
emprego potencial e o emprego efectivo no ano de 
1992, para os diferentes níveis profissionais. Verifi- 
ca-se que as carências mais acentuadas se notam nas 
categorias de empregados altamente qualificados e 
para trabalhadores manuais altamente e semi-quali- 
ficados. 

O estudo referido do ME/BM calculou os «out- 
puts» do ensino politécnico para 1992 em duas situa- 
ções: mantendo-se a política educativa actual e modi- 
ficando-se esta de modo a satisfazer as necessidades 
referidas em 1992 de acordo com os pressupostos 
apresentados (Quadro II). 

Para o ensino secundário profissionalizante que 
urge generalizar no país, temos os valores do Qua- 
dro III. 

Trata-se de valores médios de diplomados com 
os cursos profissionalizantes a criar ao nível do 9.º ano 
e com trés anos de educação média. 

Quanto ao ensino superior universitário, o rela- 
tório do Banco Mundial que temos vindo a seguir, 
diz que as taxas actuais permitirão a satisfação quan- 
titativa das necessidades de mão-de-obra durante a 
década de 80 e princípios de 90, salvo os casos que 
se referem a seguir. Assim, se não houver melhor 
redistribuição e utilização dos engenheiros agróno- 
mos, haverá excedentes no final da década de 80. 
Poderá também haver excedentes em engenharia quí- 
mica, electrónica e ciências naturais. O maior exce- 
dente, todavia, verificar-se-á entre os médicos. 

Os maiores défices serão nos seguintes casos: 


— professores do ensino secundá- 

rio e universitário — 15 a 36000 
— gestores de empresas — 12 a 16000 
— engenharia civil e arquitectura— 2,5 a 4500 


2— A LD&D NACIONAL E O DESENVOLVI- 
MENTO 


Vive-se na época da investigação científica. As 
nações que ocupam a vanguarda do desenvolvimento 
e que pontificam politicamente procuram ter e «fazer» 


cientistas. É assim que o Homo investigans é um 
recém-chegado. É uma sub-espécie do Homo sapiens 
sapiens gerado pela revolução industrial e que tem 
como ambiente privilegiado a civilização tecnotrónica 
que estamos vivendo. Não esqueçamos que estão 
vivos 90% de todos os investigadores que houve no 
mundo desde a antiga Grécia. 

O conhecimento científico e técnico em um país 
e a maneira como é usado nos processos de produção 
constituem elementos de poder que se reflectem nas 
relações desse país com os restantes. A superioridade 
ou a dependência científico-tecnológica têm signifi- 
cado político e são factores fundamentais do exer- 
cício do poder. 

A capacidade de uma sociedade em assimilar, 
adaptar, melhorar e criar conhecimento em activi- 
dades produtivas baseia-se em três pressupostos inter- 
dependentes: 

O primeiro diz respeito ao conhecimento cientí- 
fico e técnico que essa sociedade possui e ao modo 
como o usa para fins produtivos. O segundo refere-se 
à disponibilidade em recursos humanos treinados 
capazes de compreender e usar eficazmente aqueles 
conhecimentos nas suas actividades. O terceiro en- 
volve a estrutura e a eficiência das instituições ligadas 
ao progresso científico e técnico e às suas aplicações 
nas actividades económicas e sociais, nomeadamente 
as Universidades, os Laboratórios Nacionais de 
investigação científica e tecnológica, os departamen- 
tos governamentais reguladores, orientadores e finan- 
ciadores das actividades de ILD&D e as empresas. 

A aplicação desta principiologia de modo a que 
a Ciência e a Tecnologia se tornem instrumentos de 
desenvolvimento económico e social da nação, supõe 
que, o país tenha infraestrutura científico-tecnoló- 
gica adequada, suporte de capacidade científico-tecno- 


QUADRO III 


<Outputs» para satisfazer as 
necessidades de mão-de-obra 
em 1992 
(valores médios anuais) 


CURSOS 


Contabilidade e Administra- 


ção 12 000 — 15 000 
Administração Pública | 9 000 — 15 000 
Construção Civil 1700 
Química 1 600 
Máquinas 2 800 
Electricidade 1800 
Design industrial e Comercial 600 
Artes decorativas e marcena- 

ria 450 
Artes gráficas 250 
Línguas, hotelaria, turismo, 

etc. 850 
Pesca 500 


Fonte: MEU (1981) Bol. Inf. n.º 14. 
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QUADRO IV 


Evolução do número de prémios Nobel 


FÍSICA 


PAIS - 
Ante 1940 Pós 1940 
| França 6 2 
Inglaterra 10 10 
Alemanha 10 3 
Estados Unidos 6 40 
União Soviética [0 7 


Fonte: Le Monde Aujourd'hui, 14-10-1984. 


lógica. Esta capacidade exige recursos de vário tipo, 
a saber: humanos, financeiros e Institucionais. Se os 
dois primeiros tipos de recursos são fáceis de avaliar, 
os institucionais englobam os níveis do planeamento, 
da coordenação, da promoção e da realização de 
actividades de I,D&D. 

Este facto foi muito recentemente posto em evi- 
dência em um suplemento científico do jornal francês 
«Le Monde» que, para ilustrar a «fraqueza» da com- 
petitividade francesa no mundo tecno-científico publi- 
cava o Quadro IV que mostra bem como no após 
guerra a Europa tradicionalmente mãe-científica se 
deixou marginalizar. 

É que a dependência científico-tecnológica paga-se 
caro pois se traduz na perda do controlo na decisão 
da programação, produção.e mercado, na importação 
não apropriada de conhecimentos e na perda de poder 
de negociação na aquisição de tecnologia. 

Os três pontos explicitados: conhecimento cientí- 
fico e técnico disponível, recursos humanos treinados, 
estrutura institucional são pressupostos básicos que 
nos permitirão procurar perspectivar e prospectivar 
o caso português. 

A consideração dos três pontos referidos tem de 
se relacionar com a política tecnológica industrial 
e energética definida até ao fim do presente século. 
Os principais vectores de tal política aparecem men- 
cionados em documentos oficiais como o Plano de 
Desenvolvimento Tecnológico da Indústria Transfor- 
madora portuguesa e o Plano Energético Nacional. 
São eles: 


— reforço da infraestrutura científica e tecnoló- 
gica industrialmente orientada; 

— apoio à modernização, aumento de produtivi- 
dade, melhoria da qualidade e inovação dos 
sectores industriais tradicionais bem como o 
desenvolvimento de novas tecnologias para o 
uso mais racional da energia; 

— aquisição, demonstração e desenvolvimento de 
tecnologias avançadas; 

— maximização do uso dos recursos naturais na- 
cionais; 
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QUÍMICA MEDICINA 
Ante 1940 Pós 1940 | Aneis — Pós 1940 
O 3 4 
15 7 14 
ê 8 5 
23 6 51 
l 2 O 


— acompanhamento da evolução internacional das 
diversas tecnologias e suas aplicações; 

— criação de mecanismos financeiros, fiscais e 
institucionais de apoio à inovação industrial e 
uma melhor ligação entre as Instituições de 
Investigação e a Indústria. 


O Plano de Desenvolvimento Tecnológico da 
Indústria Transformadora portuguesa está dividido 
em seis programas de que apenas citaremos o pri- 
meiro: Reforço das capacidades científicas e tecno- 
lógicas com ênfase na Investigação e Desenvolvi- 
mento do sector produtivo, 

Embora este programa exija participação intensa 
do Estado, são esperadas acções relevantes de I,D&D 
nas empresas privadas. Todavia este programa supõe 
a atribuição de uma percentagem mais significativa 
do PNB às actividades de I,D&D orientadas para o 
desenvolvimento económico procurando-se atingir 
valores aceitáveis a nível dos padrões internacionais. 
Prevê-se que os gastos de investigação e desenvolvi- 
mento industrialmente orientados passem do valor 
actual de 0,12 do PNB para 0,40% no espaço de uma 
década. 

O financiamento global dos programas referidos 
terá uma componente de 60% proveniente do Estado 
e de 40% proveniente da Indústria. 

Servindo-nos dos valores globais disponíveis sobre 
recursos humanos em I,D&D publicados pela JNICT 
em 1983, verificamos que em 1980, em termos de 
equivalência a tempo integral (ETI), aqueles recursos 
tinham as seguintes distribuições (Quadro V); 


QUADRO V 


Sector de execução 


| Estado 
Ensino superior 


I.P.s/F.L. 
Empresas 


TOTAL 


LPs/F.L. = Instituições Privadas s/fins lucrativos. 
Fonte: INICT (1983) Est. Doc. Ciên. & Tecnologia n.º 2. 


Se agora procurarmos discriminar estes 7710,9 ETI 
por funções exercidas, temos os valores do Qua- 
dro VI. 

Verifica-se que 34,5% correspondem à função de 
investigador, 37,2% são técnicos e os restantes 28,3% 
são operários. 

Um facto relevante a reter da análise dos Qua- 
dros V e VI refere-se à modéstia, em termos rela- 
tivos, dos efectivos de LD&D no sector Empresas 
devido ao discreto número de investigadores quando 
comparados com os números correspondentes dos 
sectores do Estado e Ensino Superior. 

Averiguando quais as qualificações dos recursos 
humanos em LD&D obtêm-se os valores do Qua- 
dro VII 

Este quadro permite ter uma visão bimodal dos 
recursos humanos adstritos à LD&D. Assim temos a 
distribuição percentual pelos quatro sectores aí refe- 


ridos e ainda a distribuição por qualificações. 

Os três quadros apresentados (Quadros V, VI e 
VII) podem sintetizar-se em termos percentuais como 
se aponta no Quadro VIII. A complementá-lo estão os 
valores do Quadro IX em que os recursos humanos 
aparecem desagregados por domínios científicos e 
tecnológicos. 


Valerá a pena considerar o Quadro X onde se 
fazem algumas comparações entre países europeus. 

Para se ter uma noção mais exacta das assimetrias 
da I,D&D no nosso país convirá atentar nos valores 
do Quadro XI onde o país aparece dividido em qua- 
tro grandes áreas. 


Procurando mostrar um panorama razoavelmente 
actual, geral e funcional de recursos humanos nacio- 
nais apresentam-se os quadros seguintes (Quadros XII 
a XV). Aí aparecem a distribuição dos recursos finan- 


QUADRO VI 


| Sector de | 
execução | Estado 
Função | 


Ensino Superior 


Empresas 


Investigadores 836,2 1 380,2 


Técnicos 1 678,4 549,0 
| Operários | T21,8 289,1 
Administrativos 474,9 160,0 
TOTAL STLio 23783 
Fonte: idem Quadro V. 
QUADRO VII 


= o Sector de 
execução 
Qualificação | 


Curso Superior 
| Curso médio 
| Outros * 


TOTAL 


48,1 | 
“* Vidé nota Quadro VIII. 


Fonte: idem Quadro V. 


QUADRO VIII 


Ensino Superior I.P.s/F.L. 


Estado Empresas 
| o io % 

FUNÇÃO 

Investigadores | 22,5 580 35,2 260 

Técnicos 45,2 23,1 320 | 41,0 

Operários | 19,5 12,2 15,3 

Administrativos | 12,8 6,7 17,5 a 
QUALIFICAÇÃO | 

Curso superior 39,2 65,2 60,3 | 26,4 

Curso médio 28,6 16,0 | 69 | 40,5 

Outros * 322 


* Pessoal com estudos secundários incompletos, operários especializa- 


dos ou não. 


18,8 | 32,8 33,1 


Fonte: idem Quadro V. 
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QUADRO IX 


Domínio Científico 
RA] 
Tecnológico 

| | 
Ciên. Exactas e Naturais 
Ciên. da Engenharia 
Ciên. da Saúde 
Ciên. da Agricultura, Silvicultura, 
Pecuária e Pescas 
Ciênc. Sociais e Humanas 


Fonte: idem Quadro V. 


Ensino Superior 


Empresas 


So 


% 


QUADRO X 
Pessoal de I&D em ETI, dados referentes a 1977 


Pessoal em I&D (ETI) 


Total 


Graduados Universitários 
“% do total 


0/o da força laboral 


R. Fed. da Alemanha 319 434 
Bélgica 29 532 
Dinamarca 13 801 
França 222 111 
Holanda 52 857 
Irlanda 5 857 
| Itália 97 345 
Portugal 6 543 


12 23 
7 31 
5 39 
10 não disponível 
1 28 
5 não disponível 
5 37 
2 44 


Fonte: Situação do sistema científico e tecnológico português MCCC (1982). 


QUADRO XI 
I, D&D por regiões do país (1978) 


Regiões Pessoal (ETI) 


bri 
Distritos Entre-Douro e Minho 
Distritos Entre-Douro e Tejo 
(salvo Lisboa) 
Grande Lisboa 
Alentejo, Algarve e regiões autó 
| nomas dos Açores e Madeira 


Fonte: KIM, H.K. (1982). 


ceiros para 1983 que contemplaram os principais 
laboratórios nacionais estatais (Quadro XII) e a ilus- 
tração dos recursos humanos por ramos de actividade 
(Quadros XIII a XV). Os valores apresentados são 
eloquentes e falam por si, mais do que os glosar 
interessa que meditemos neles. 


3 — CONSIDERAÇÕES FINAIS 


É importante, após a quantificação que se pro- 
curou ao juntar os quadros apresentados, referir que 
o Plano Tecnológico identificou as áreas programá- 
ticas de alta prioridade que representam as necessi- 
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dades e oportunidades tecnológicas para a indústria 
portuguesa. 

São elas: 

— Necessidade de aumentar e fortalecer a capa- 
cidade analítica nacional necessária ao desen- 
volvimento de uma política tecnológica. 

— Necessidade generalizada de melhorar a capa- 
cidade tecnológica das empresas existentes. 


QUADRO XII 


Recursos financeiros dos organismo nacionais 
de I, D&D em 1983 (10º escudos) 


Ordinário 


LNETI (Energia e Tecnologia) 
LNEC (Engenharia Civil) 
INIAER (Agricultura) 
INIC (Inv. universitária) 
HCT (Inv. tropical) 
INIP (Pescas) 
IH (Hidrografia) 
LNIV (Inv. veterinária) 
| INICT (Política científica 
e tecnológica) 


| TOTAL 


Fonte: Examens des politiques nationales de la science. 
Portugal. OCDE (1984). 


— Necessidade de reformulação da estrutura fis- cidades nacionais nos domínios da metrologia, 


cal mais apropriada à incentivação do desen- microelectrónica, ciências dos materiais e na 
volvimento tecnológico. biotecnologia. 
— Necessidade de introduzir novas tecnologias e — Necessidade de melhorar o capital humano que 


novas firmas de base tecnológica. está na base do processo de desenvolvimento 


— Necessidade em apostar no aumento das capa- tecnológico. 


QUADRO XIII 
Recursos humanos em I, D&D, por função e por ramo de actividade económica (1982) 


Técnicos (!) 
ETI 


Investigadores | 


RAMO DE ACTIVIDADE ETI 


Extracção de minérios metálicos 
Alimentação, bebidas e tabaco 

Madeira e Cortiça 

Papel e Artes Gráficas 

Ind. Químicas dos derivados do petróleo e do carvão 
Produtos farmacêuticos 

Ind. prod. minerais não metálicos 

Ind. Metalúrgicas de base 

Fab. produtos metálicos 

Máquinas não eléctricas 

Máquinas e aparelhos eléctricos 

Material de transportes 

Aparelhos de medida, fotográficos e óptica 
Electricidade, gás e vapor 

Construções e Obras Públicas 
Comunicações 

Serviços prestados às Empresas 


(1) Com curso secundário ou equivalente. 


Fonte: I&DE. Dados de actualização a 31-12-82. Valores provisórios. JNICT (1984). 


QUADRO XIV 


Recursos humanos em actividades de Controlo de Qualidade, por função 
e por ramo de actividade económica (1982) 


Técnicos com Técnicos 


RAMO DE ACTIVIDADE Curso LN ETI ETI 
Ext. de minérios metálicos 20 15,0 
Alimentação, bebidas e tabaco 67,06 122,2 226,64 
Madeira c Cortiça 22,1 180,9 217,5 
Papel e Artes Gráficas 4,0 7,67 
Ind. Quim. deriv. petról. e carvão 177,61 515,89 
Produtos farmacêuticos 40,0 91,0 
Ind. prod. min. não metálicos 80,47 102,57 
Ind. Metalúrgicas de base 369,0 380,0 
Fab. produtos metálicos -— 5,0 
Máquinas não eléctricas 31,5 62,11 
Máquinas e aparelhos eléctricos 130,5 264,75 
Material de transportes 0,5 84,5 
Aparelhos medida, fotog. e óptica 35,0 45,98 
Electricidade, gás e vapor 24,0 198,0 
Construções e Obras Públicas — 6,0 
Comunicações 30 223 
Serviços prestados às Empresas 2,0 76,9 


1 235,13 2350,66 | 


322,72 


Fonte: I&DE. Dados de actualização a 31-12-82. Valores provisórios. 
INICT (1984). 
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QUADRO XV 


Recursos humanos em actividades científicas excluindo Controlo de Qualidade, 
por função e por ramo de actividade económica (1982) 


RAMO DE ACTIVIDADE 


Extracção de minér. metálicos 
Alimentação, bebidas e tabaco 
Madeira e Cortiça 
Papel e Artes Gráficas 
Ind. Quím. deriv. petróleo e carvão 

| Produtos farmacêuticos 

| Ind. prod. min. não metálicos 

| Ind. Metalúrg. de base 

' Fab. produtos metálicos 

' Máquinas não eléctricas 
Máquinas e aparelhos eléctricos 
Material de transportes 
Apar. med., fotog. e óptica 
Electricidade, gás e vapor 
Const. e Obras Públicas 
Comunicações 
Serviços prestados às Empresas 


TOTAL 


| Curso superior 


Técnicos com dis 
Auxiliares 


ETI é 


12,16 
19,31 
1,9 
34 
2,65 
131,65 
8,0 
0,64 
35,76 
3,05 
3,75 
37,0 
7,0 
83,0 
2,0 
92,25 
257,3 


700,82 2 326,11 


Fonte: I&DE. Dados de actualização a 31-12-82. Valores provisórios. 


INICT (1984). 


Devem estas seis necessidades básicas ser satisfei- 
tas para que o país entre na via de um desenvolvi- 
mento harmonioso. O suporte humano é essencial, 
para obviar a estas necessidades e há que aproveitar 
as oportunidades que elas nos põem, procurando: 


— Criar um conjunto de organizações que acom- 
panhem o desenvolvimento tecnológico em Por- 
tugal e no estrangeiro, identificando as neces- 
sidades e oportunidades tecnológicas e desen- 
volvendo as acções oportunas. 

— Melhorar a produtividade e a qualidade dos 
produtos e o desenvolvimento de novas tecno- 
logias através da actividade de normalização e 
certificação de qualidade, do estabelecimento de 
centros tecnológicos e de uma melhor colabo- 
ração entre as universidades, os laboratórios 
nacionais e a indústria promovendo uma crite- 
riosa transferência de tecnologia. 

— Promover a inovação inclusivamente pela cria- 
ção de um organismo, que a incentive, tipo 
agência para o desenvolvimento de tecnologias 
de ponta. 

— Difundir a informação técnica de modo a pro- 
vocar o arranque de empresas de base tecno- 
lógica acentuada. 


— Perspectivar o apoio no uso dos recursos natu- 
rais do país procurando cobrir campos cruciais 
para o nosso desenvolvimento no contexto da 
inserção no mundo fortemente tecnotrónico. 


— Formar pessoal técnico a todos os níveis, ges- 
tores de tecnologia e trabalhadores especializa- 
dos, tudo gente capaz de responder aos desa- 
fios da civilização não só actual, mas que se 
antevê para amanhã. 


O esforço para melhorar a qualidade dos recursos 
humanos do país terá que ser muito grande. Por um 
lado, precisamos de empresários e gestores profissio- 
nais com formação em engenharia ou pelo menos 
técnica, capazes de gerir actividades inovadoras. Por 
outro lado urge incrementar a preparação dos traba- 
lhadores nos seus locais de trabalho de modo a 
aumentarmos o número de técnicos e de operários 
especializados. Só procurando a satisfação dos pontos 
acima mencionados, a transferência de tecnologia será 
via transitória para independência nacional e não via 
permanente para subordinação ao exterior. 


Recordamos que uma das causas determinantes 
da nossa falta de desenvolvimento tecnológico está na 
escassez de I,D&D a suportá-lo. Acresce que só cerca 
de 13% da I,D&D portuguesa é de iniciativa do sector 
industrial privado enquanto que cerca de 66% é de 
iniciativa do governo. Mostram os factos que nos 
países industrializados a I,D&D privada está mais 
ralacionada com o crescimento económico que a 
I,D&D estatal. 


Um estudo envolvendo cem empresas nacionais 
de diferentes tamanhos e referente a 1981 e 1982 
mostra que estas empresas se repartiam pelos seguin- 


tes níveis tecnológicos de acordo com os valores per- 
centuais do Quadro XVI. 


QUADRO XVI 


Percentagem 


Fonte: Plano de Desenvolvimento Tecnológico da Indústria 
Transformadora Portugesa. Conclusões e Recomen- 
dações LNETI (1983). 


Os níveis são caracterizados como segue: 


Nível |: Empresas que têm actividades de I&D e 
engenharia (com departamentos próprios), 
que podem desenvolver novos produtos € 
processos e introduzir alterações importan- 
tes em produtos e processos já existentes. 

Nível 2: Empresas só com departamentos de enge- 
nharia e que podem introduzir modifica- 
ções importantes em produtos e processos 
(geralmente com tecnologia adquirida no 
exterior). 

Nível 3: Empresas com pouco pessoal técnico e que 
podem apenas introduzir pequenas altera- 
ções em produtos e processos (com tecno- 
logia adquirida). 

Nível 4: Empresas com orientação tecnológica mas 
sem capacidade interna para além de ope- 
rarem competentemente a tecnologia que 
possuem (se uma fábrica subsidiária de uma 
multinacional for indexada ao nível 4 isto 
pode querer dizer, por exemplo, que apenas 
procede à montagem dos seus produtos em 
Portugal, com tecnologia importada, de 
forma a tirar partido do preço vantajoso 
de mão-de-obra). 


Nível 5: Empresas nos sectores tradicionais e sem 
grande interesse por melhoramentos tecno- 
lógicos. 


A vantagem que a economia nacional tem tido, 
decorrente da mão-de-obra barata, tende a esbater-se. 
Resulta daqui que a rendibilidade da produção será 
o factor determinante da competitividade dos novos 
produtos nos mercados internacionais. A desvantagem 
em relação aos países estrangeiros só pode ser ven- 
cida com o aumento da capacidade inovadora do país. 

Entende-se por inovação a aplicação de uma inven- 
ção ao processo produtivo ou à estruturação de um 
serviço. Ou dito de outra forma, é a arte de associar 


elementos pré-existentes para criar novas formas, 


produtos ou processos comercializáveis. A inovação, 
exige enorme esforço de gente capaz de tomar as 
invenções e as passar ao processo produtivo e em tal 
acção promover o desenvolvimento da nação. 

Entende-se por desenvolvimento a evolução glo- 
bal de uma colectividade — no geral uma nação — 
para uma direcção positiva permitindo a satisfação 
equitativa das necessidades físicas, psíquicas, culturais 
e morais das gentes. Esta perspectiva de desenvolvi- 
mento supõe que o Homem é o valor supremo e 
o fim da História e que a Ciência e a Técnica têm 
funções meramente instrumentais. 

Nesta perspectiva o fortalecimento dos recursos 
humanos do país é a prioridade das prioridades. 
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DIAGRAMAS POTENCIAL — PH 
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RESUMO 


Uma das contribuições mais importantes das uúlti- 
mas décadas para a compreensão dos fenómenos de 
corrosão metálica deve-se a Pourbaix e seus colabo- 
radores que construiram diagramas de equilíbrio ter- 
modinâmico (potencial-pH) para o estudo sistemático 
aa previsão das reucções de metais em soluções aquo- 
sas. Neste artigo estabelecem-se estes diagramas por 
uma via simplificada, baseada no uso directo da 
equação de Nernst, método que se tem mostrado 
particularmente adequado na introdução do assunto 
aos estudantes. Faz-se ainda referência ao emprego 
e limitações dos diagramas, bem como a alguns exem- 
plos de sistemas de carácter industrial. 


ABSTRACT 


One of the most important contributions to corro- 
sion science and engineering in recent decades has 
been the work of Pourbaix and associates in the deve- 
lopment of termodynamic equilibrium diagrams (po- 
tential vs pH) for the behaviour of metals in aqueous 
solutions. The author presents herein a simplified 
procedure for constructing these diagrams, based on 
the direct use of the Nernst equation, which has 
been successful in presenting the subject to students. 
Applications and limitations of the Pourbaix dia- 
grams, as well as some examples of systems with 
technological interest are also given herein. 


Original recebido para publicação em 22/3/84. 
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INTRODUÇÃO 


Um dos problemas mais importantes que se põem 
em Corrosão, consiste na selecção de materiais téc- 
nicos que vão ser usados em contacto com meios 
aquosos, em processos tecnologicamente importantes. 
A aprovação ou recusa de um determinado material] 
deve decidir-se em função da velocidade de dissolução 
desse material no meio em que vai ser usado. Em 
geral, porém, e para evitar tentativas inúteis e anti- 
-económicas, é costume começar por se recorrer à 
termodinâmica pois esta permite prever o sentido 
possível das reacções de dissolução e de redução, que 
constituem o processo corrosivo. Note-se que o mé- 
todo termodinâmico não permite senão prever, isto é, 
afirmar quando uma reacção é possível ou impossível, 
sem poder garantir, naquele caso, que ela realmente 
se dê. Contudo estes conhecimentos são, por si só, 
suficientemente importantes e úteis para a abordagem 
inicial do problema da selecção de materiais. 

O estabelecimento do sentido da evolução provável 
de uma reacção electroquímica, por via termodiná- 
mica, pode atingir-se por diversos caminhos, mas é 
vulgar escolher como grandezas que melhor se pres- 
tam para representar as condições de evolução das 
reacções electroquímicas e as respectivas funções de 
estado, o potencial (ou tensão) reversível de elécrodo 
e a energia livre duma reacção. Fixa-se então o sen- 
tido da evolução daquelas pelo sentido de variação 
da energia livre que as caracteriza. Se AG < O para 


um processo corrosivo, então a corrosão do metal 
verifica-se espontaneamente. Isto equivale a dizer que 
a tensão (ou potencial) do sistema metal/meio subiu 
acima do seu valor de equilíbrio. Efectivamente, se 
G->G', tal que AG=G'—- G=nF (E- EJ<O0, 
então E>E' cnE=E -—-E>O, cq.d. 

Para os sistemas electroquímicos em que um ou 
mais dos seus constituintes estão envolvidos em diver- 
sos tipos de equilíbrio, é extremamente útil recorrer 
aos métodos gráficos porque é então possível estabe- 
lecer o diagrama referente a todas as condições pro- 
váveis de evolução dos constituintes do sistema 
perante todo o conjunto de circunstâncias que afec- 
tam os factores internos dos equilíbrios, desde que os 
factores externos (P e T) se mantenham constantes. 
Diversos diagramas têm sido estabelecidos por diver- 
sos modos. Dentre estes destacam-se os diagramas 
devidos a Marcel Pourbaix [1] que se destinam essen- 
cialmente a pôr em evidência a influência dos iões H* 
sobre os equilíbrios. Inúmeros são os problemas que 
os diagramas permitem resolver, de Química Analí- 
tica, Electroquímica, Química Inorgânica, etc., mas 
a sua principal importância resulta da forte contri- 
buição que trouxeram ao estudo da Corrosão Metá- 
lica. Neste domínio, como é óbvio, o estudo da 
influência dos iões H* sobre os equilíbrios reveste-se 
de especial significado prático visto que grande parte 
dos sistemas electroquímicos de interesse tecnológico 
são do tipo metal/solução aquosa, em que estão pre- 
sentes os iões H' os quais, dadas as propriedades que 
lhes são peculiares, têm acção preponderante sobre 
a tensão de eléctrodo. As representações gráficas de 
Pourbaix baseiam-se no princípio termodinâmico 
exposto inicialmente; tais representações permitem 
englobar sistematicamente num único diagrama, todos 
os estados de oxidação-redução referentes a todas as 
formas oxidadas e reduzidas, que um mesmo elemento 
químico pode formar com dado reagente. Estes dia- 
gramas são construídos de modo que cada linha —- 
uma recta — represente condições de equilíbrio entre 
duas formas, apenas, do sistema global; estes equilí- 
brios podem ser do tipo ácido-base ou do tipo oxida- 
ção-redução. Resulta que o encontro de linhas com 
formas comuns, delimita zonas de estabilidade de 
cada uma delas; esta constitui uma das particulari- 
dades mais importantes destes diagramas. 


REPRESENTAÇÕES GRÁFICAS 


A construção dum diagrama E-pH, a 25ºC, im. 
plica a representação de equações do tipo: 


oro DOI O 005 | 
== o + — es TT 
nº YR E 
Cc 
log — = log K —b pH 
E g 


que caracterizam, respectivamente, equilíbrios elec- 
troquímicos do tipo: 


a0+bHt+ne=cR+dHO 
e equilíbrios químicos do tipo: 
a A+b Ht=c D+ dHO 


onde b e d são geralmente inteiros e dígitos podendo 
ser positivos ou negativos; a, c e n são sempre núme- 
ros positivos geralmente inteiros e quase sempre 
dígitos. Note-se que se fizeram intervir os electrões 
no 1.º membro da reacção electroquímica a fim de 
respeitar a Convenção Europeia de Sinal, convenção 
essa que é adoptada ao longo deste artigo. Observe-se 
ainda que se considerou a actividade da água cons- 
tante, isto é, supõe-se tratar-se de soluções suficiente- 
mente diluídas. 

A representação linear das equações escritas acima 
faz-se fixando previamente uma relação, arbitraria- 
mente escolhida, para os quocientes (0) /(R)J e 
(DJ / (A). Há portanto famílias de rectas correspon- 
dentes aos diversos valores destes quocientes. Esta 
representação oferece a vantagem de poder traduzir 
as variações de E (ou de pH) quando se fazem variar 
aqueles quocientes e ainda a de mostrar que cada 
zona do plano E-pH delimitada por rectas represen- 
tativas de equilíbrios em que intervém um mesmo 
constituinte representa também além da estabilidade 
a predominância desse constituinte. Cada ponto do 
plano possui três coordenadas, duas de referência aos 
eixos ortogonais e uma terceira em relação às rectas 
fronteiras de cada domínio; assim uma forma predo- 
minará tanto mais sobre as outras quanto mais afas- 
tado se encontrar o ponto representativo do sistema 
das linhas fronteiras do domínio. Com efeito, consi- 
derando a concentração 10º moles/litro como limite 
de detecção analítica corrente para as formas dissol- 
vidas, a distância entre cada linha correspondente à 
variação de 10 ursdades das concentrações molares, 
é de 0,0591/z volt, para os 10es de valência z à tem- 
peratura de 25º C, donde resulta que a variação das 
actividades desde | molar até 10% molar ocupe dimen- 


a 


sões do plano geralmente restritas em relação à área 
de cada domínio. A zona correspondente à estabili- 
dade dum dado constituinte delimitada pelas linhas 
de concentração 10% molar, representa pois também 
c domínio de predominância da forma respectiva. 

Feitas estas considerações representemos sucessi- 
vamente no plano E-pH as linhas correspondentes aos 
diversos equilíbrios susceptíveis de ocorrer nos siste- 
mas do tipo metal-solução aquosa. 


Reacções electroquímicas dependentes do pH 


A representação das rectas correspondentes aos 
equilíbrios de oxidação-redução em que intervêm os 
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iões Ht ou OH" pode fazer-se por uma família de 
linhas oblíquas paralelas cuja inclinação será função 
do coeficiente estequiométrico que afecta o H* ou 
o OH” na reacção e do número de electrões postos 
em jogo, e cuja posição face ao eixo dos potenciais 
dependerá da actividade dos diferentes constituintes, 
ionisados, ou não, que tomam parte na reacção. 

Como exemplos destas reacções podemos apontar 
os seguintes: 


MnoO; + 8H* + Se = Mn” + 4H,0 
Fe(OH), + 3H+ + e = Fe” + 3H,0 
2H + 2e = H, 
(MnO)) 
= 7— € OO; esse—m=iio 
E = + 1,507 — 0,0946 pH + 0,0118 log (Mn) 
E = + 1,057 — 0,1773 pH— 0,0591 log (Fe*) 
E = 0,000 — 0,0591 pH — 0,0295 log py, 
Observe-se que estas famílias de linhas se reduzem 
a linhas apenas nos casos em que nas reacções não 


participam constituintes gasosos nem dissolvidos para 
além dos iões H*, como acontece com o equilíbrio: 


Fe(OH), + 2H* + 2e = Fe + 2H,0 
E = — 0,047 — 0,0591 pH 


EV 


o Mn; 
e 
[MO] 
q [Mind] 
a (Mnaç 
Mn! | [Mint] 
Fig 1 Representação de equilibrios 
eleciroquimicos dependentes 
pm do ph. 


Na Fig. 1 representa-se a família de linhas corres- 
pondentes ao equilíbrio MnO,/Mn*. Estas linhas divi- 
dem o plano E-pH em duas regiões: a de estabilidade 
termodinâmica dos iões MnQO, e a de estabilidade ter- 
modinâmica dos iões Mn?, Note-se que se os iões H* 
aparecessem no membro da reacção onde se encontra 
a forma reduzida (Mn”), o poder oxidante do sistema 
aumentaria com o pH. 


Reacções electroquímicas independentes do pH 


A representação das rectas correspondentes aos 
equilíbrios electroquímicos independentes do pH pode 
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fazer-se por uma família de linhas horizontais cuja 
posição dependerá da actividade dos diferentes cons- 
tituintes, ionisados ou não, que tomam parte na reac- 
ção. Como exemplo destas reacções podemos apontar 
as seguintes: 


Fer + é = Pe” 


Fe“ + 2 e= Fe 
CL+2e=2 CI 


(Fe!) 
E = + 0,771 + 0,0591 log —— 
E “(rem) 
E = — 0,440 + 0,0295 log (Fe”) 
Pci, 
E = + 1,359 + 0,0295 log E 
Cl 


= fio 
«Qi e mma 
- jrero 
Fe 
Fig 2 Representação de equilibrios 
mio plectroquimicos independentes 
pH do pH. 


Na Fig. 2 representa-se a família de linhas corres- 
pondente ao equilíbrio Fe?*/Fe. Estas linhas dividem 
o plano E-pH em duas regiões: a de estabilidade ter- 
modinâmica dos iões Fe” e a de estabilidade termo- 
dinâmica do sólido Fe. 


Reacções químicas dependentes do pH 


A representação das rectas correspondentes aos 
equilíbrios químicos que dependem do pH pode fazer- 
-se por uma família de linhas verticais cuja posição 
dependerá da actividade dos constituintes, ionisados 
ou não, que tomam parte na reacção. 

Como exemplo destas reacções podemos apontar os 
seguintes: 


HCO; + H* = HCO, 
Fe(OH), + 2 Ht = Fe” + 2H,0 
HCO; + H* = CO, + HO 


lo ARGOS =—6,37+pH 
é (HCO; RES 
log (Fe*) = 13,29 — 2pH 
Pco, 
log ———— = 77,80 — pH 
(HCO,) 


